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Chapter 1 



Introduction 



Since the discovery of the neutron, immediately fohowed by their prediction, neutron 
stars are one of the more fascinating and exotic bodies in the universe. Matter inside 
reaches densities very far to be reahzed in present experiments, and theoretical 
prediction can be verified only indirectly with astronomical observations. 

Neutron stars [U [H [3l H] belong to a particular class of stars whose mass is 
generally of the order of ~1.5 solar masses (Mq) embodied in a radius of ~12 Km, 
and have a central density ranging from 7 to 10 times the nuclear saturation density 
po=0.16 fm~^ which is found in the center of heavy nuclei. The Fermi energy of 
Fermions at such densities is in excess of tens of MeV and then thermal effects are a 
very small perturbation to the structure of neutron stars. Therefore they exhibit the 
properties of cold matter at extremely high densities and have proved to be fantastic 
test bodies for theories from microscopic nuclear physics to general relativity [5]. 

The crust of neutron stars (the outer 1-2 Km) primarily contains nuclei in the 
region where the density increases toward the crust-core interface at p ~ Po/S, 
and is made in large part of nuclei in which the proton fraction is x ~ 0.1 to 0.2; 
such extremely neutron-rich nuclei cannot be directly observed, but can be created 
in laboratories in heavy-ion beam experiments [6j. At higher densities neutron rich 
nuclei undergo a number of transitions through phases known as nuclear pastafUPT], 
until the homogeneous neutron matter phase is reached. 

When the mass of neutron-rich nuclei exceeds the so called drip-line, neutrons 
form a uniform gas. At the crust-core interface the mixing of nuclei and neutron 
matter gives rise to a phase in which nuclei are probably organized on a lattice. 

The neutron matter in the inner crust may be in a ^5*0 superfluid phase[8] that 
alters the specific heat and the neutrino emission of the star playing then a crucial 
role in the star cooling. The superfluid phase is also investigated in order to possi- 
bly explain sudden spin jumps observed in neutron stars, called glitches, probably 
related to the angular momentum stored in the rotating superfluid neutrons in the 
inner crust [1]. 

The outer core of the neutron star consists of matter composed by nucleons. 
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electrons and muons. Nucleons are in the larger part neutrons with a small fraction of 
protons, and the degenerate electron gas contrasts the gravitational field by avoiding 
that the star becomes a black hole. 

In the inner core higher hadrons could appear and play a crucial role in the 
matter (called then hadronic matter). Hadrons could be charged mesons, such pions 
or kaons that eventually could form a Bose condensate[9l [10], and other hyperons 
such as the S~ or A. The nucleon-hvp er on | IH [T2] and hyp er on-hyp er on } 1 3j inter- 
actions are not as well known[14] as those between nucleons. As a consequence the 
transition from nucleonic to hadronic matter is hard to be calculated. Some at- 
tempt to explicitly include hyperons in the nuclear matter equation of state was 
recently proposed [El [16]. At higher densities there should also be a transition to 
quark matter in which quarks are deconfined |17j . 
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The simpler model to predict observable interesting properties of neutron star 
starts with the study of pure neutron matter bulk, and after with the addition of 
a small fraction of protons. The modelization and the knowledge of the nucleon- 
nucleon interaction from scattering data is then the starting point for a realistic 
simulation of neutron stars. At the same time the research of an accurate technique 
to exploit such effective interactions is fundamental to attack ab-initio problems of 
interest for nuclear astrophysics. 

This thesis is concerned about of the interesting properties of neutron matter 
briefly described above. It will be shown how the Auxiliary Field Diffusion Monte 
CarloflHlIiniEOlET] (AFDMC) method can be used to accurately study a wide range 
of phenomena that govern the physics of a neutron star. The only needed starting 
point is the choice of the degrees of freedom and consequently the Hamiltonian of 
the system. 

In nuclear physics a calculation typically starts with the assumption that nu- 
cleons are nonrelativistic particles interacting with some effective potentials. This 
approach is a great simplification over reality. Nucleons are composite systems with 
an internal structure due to quarks interacting by gluon exchange, but presumably 
the quark degrees of freedom become important only at very short distances with 
high energies and momenta exchanged, so a nucleon-nucleon potential description 
should be adequate for the realm of low-energy nuclear physics. 

At the present there are several forms of realistic Hamiltonians, able to accu- 
rately reproduce scattering processes and light nuclei properties. Unfortunately the 
development of modern potentials and of the few body-methods used to study very 
light nuclei, was not followed by the development of an accurate method to test 
and apply the nuclear Hamiltonians to the study of heavier nuclei, neutron and nu- 
clear matter, and the study of properties of four-body nuclei is the standard field of 
investigation [221 ES] . 

Given a nuclear Hamiltonian, accurate calculations are limited to very few nucle- 
ons. The exact Faddeev-Yakubovsky equation approach is limited to systems with 
A=4[24j: other few-body variational calculations [25 j solve the nuclear Schrodinger 
equation in good agreement with the exact method but are extended only up to 
A=6 nuclei [26| I27j. Techniques based on shell models calculations like the ab initio 
no-core shell model were extended up to A=13[28j and very recently up to A=40[29j 
(but with some commented uncontrolled approximation [30] ) . Quantum Monte Carlo 
methods, like the Variational Monte Carlo[31', '32] (VMC) or other techniques based 
on recasting the Schrodinger equation into a diffusion equation, like the Diffusion 
Monte Carlo (DMC) or Green's Function Monte Carlo[33l[Sl (GFMC) allowed for 
performing calculations of nuclei up to A=12[35j and neutron matter with A=14 
neutrons in a periodic box[36l [371 [38] . 

Other many-body theories that contain uncontrolled approximations are used 
for heavier nuclei. The first class of them essentially modify the Hamiltonian to an 
easy-to-solve form like the Brueckner-Goldstone|39j or the Hartree-Fock[40j methods 
and all the theories constructing effective interactions from microscopic ones. The 
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second class include methods that work on the trial wave function typically used in 
variational calculations like the Correlated Basis Function (CBF) theory [4T | 142 ^ [43] . 
Variational Monte Carlo [3H I45j or Coupled Cluster Expansion |46j. 
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Figure 1.2: Schematic representation of the phase diagram temperature-density of 
matter. pQ is the saturation density at T=0, extrapolated from the density of the core 
of massive nuclei. The diagram shows the region in which the AFDMC calculations 
were performed, for densities that range from the crust and part of the interior of 
neutron stars to ordinary nuclei at zero-temperature. 



It will be shown that AFDMC can be used to simulate both nuclei [47J and 
nuclear matter [l8] with high accuracy. After the choice of the Hamiltonian for each 
different system considered, AFDMC is applied to study properties of light nuclei 
and confined neutron drops to test the agreement of this algorithm with other exact 
methods, and then extended to heavier nuclei never studied with the same accuracy. 
This step is necessary to demonstrate that an accurate many-body calculation is now 
possible for realistic modern Hamiltonians. 

It will be then illustrated the application of AFDMC to study the symmetric 
nuclear matter equation of state, and to verify the limitations of other methods to 
deal with the more simple nucleon-nucleon interaction that contains the tensor-r 
force generated by the one pion exchange process. 

The AFDMC method was then applied to the study of the equation of state 
of neutron matter, that is the simpler starting point to study the structure of a 
neutron star. The global structure of a neutron star, given in terms of the mass- 
radius relation, is determined by the equations of hydrostatic equilibrium. For a 
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non-rotating spherical object in general relativity, these equations are called Tolman- 
Oppenheimer-Volkov|49j (TOV) equations. The numerical solution of the TOV 
equations gives the M-R relation between the total mass versus the radius, or gives 
the total mass as a function of the central density pc of the neutron star. 

By means of the AFDMC we studied the neutron matter in the lower density 
regime where neutrons form a ^5*0 superfluid phase[50]. It will be shown how the 
BCS phase of the system lower the energy with respect to the normal neutron gas, 
and therefore the superfluid energy gap is evaluated. 

In the last part of the thesis the AFDMC calculation of neutron-rich isotopes in 
a very simple model, will be presented. In this case the external neutrons interact 
with a realistic interaction but all the excitation effects of the core interacting with 
valence nucleons are neglected [51] 152] , The obtained results are in a very good 
agreement with experimental data, qualitatively better than other methods. 



Chapter 2 



Hamiltonian 



The goal of nuclear physics, and consequently of nuclear astrophysics, is to un- 
derstand how nuclear binding, stability and structure arise from the interactions 
between individual nucleons. To achieve this goal the starting point consists in the 
determination of the Hamiltonian to be used. 

In principle the nuclear Hamiltonian should be directly derived from quantum 
chromodynamics (QCD) but the direct calculation starting from QCD is at present 
too complex to be attacked. Thus the structure of a nuclear Hamiltonian is usually 
determined phenomenologically and then fitted to exactly reproduce the properties 
of few-nucleon systems[53j. However, in some case, its form is determined by the 
many-body methods used to solve the Schrodinger equation |54j. 

In this chapter the modern nuclear Hamiltonians will be illustrated, while the 
many-body techniques used to solve the nuclear ground-state will be the subject of 
the following chapters. 

Properties of a generic nuclear system can be studied starting from the non- 
relativistic Hamiltonian 

^ = -|^E^' + E%+ E + - (2.1) 

i i<j i<j<k 

which includes the kinetic energy operator, a two-nucleon (NN) interaction Vij, a 
three- nucleon (TNI) interaction Vijk, and in principles could include higher many- 
body interactions that however are expected to be less important. It has been shown 
that in light nuclei the three-body force contributes only a few percent of the total 
potential energy, while for heavier nuclei the contribution of TNI could be of about 
15-20%. One expects that in the physics of nuclei a possible four-body interaction 
should be either negligible or just a very small perturbation to the system. This is 
not necessarily true in neutron and nuclear matter in the typical density regime of 
neutron star where the contribution of many-body interactions in the Hamiltonian 
could be more important. 
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The NN interactions are usually dependent on the relative spin and isospin state 
of the nucleons and therefore written as a sum of several operators. The coefficients 
and radial functions that multiply each operators are adjusted by fitting experimen- 
tal scattering data, and the type and the number of these operators depend on the 
interaction. 

A large amount of empirical information about the nucleon-nucleon (NN) scat- 
tering problem has been accumulated over time till 1993, when the Nijmegen group 
analyzed all NN scattering data below 350 MeV published in a regular physics 
journal between 1955 and 1992 155]. NN interaction models that fit the Nijmegen 
database with a "x^ /N^ata ~1 are called "modern" and the more diffuse of them in- 
clude the Nijmegen models[56] (Nijm93, Nijm I, Nijm II and Reid-93), the Argonne 
models [571 ES] and the CD-Bonn[59]. In the last year also modern NN interaction 
derived directly from chiral effective field theory were proposed [6(2 |6T]. However all 
these NN interactions had been proved to underestimate the triton binding energy, 
by suggesting that the contribution of a TNI force is essential to reproduce the 
physics of nuclei. 

The TNI contribution is mainly attributed to the possible A intermediate states 
that an excited nucleon could assume after and before exchanging a pion with other 
nucleons. This process can be written as an effective three-nucleon interaction, and 
its parameters are fitted on light nuclei [62| [63] and eventually on properties of nuclear 
matter, such the empirical equilibrium density and the energy at saturation [64j. The 
TNI must accomplish the NN interaction and the total Hamiltonian has then to be 
studied after the choice of NN[65], then different form of TNI also in the momentum 
space were proposed [66]. 

The possibility to explicitly include the A degrees of freedom to avoid the addi- 
tion of the TNI in the Hamiltonian were explored by Wiringa et al. with the Argonne 
AV28 potential[67], but its application to calculate properties of triton revealed to 
be very difficult (results took 6 papers! [68l [69| iTO t Ffll l72l l73] ) . and also in nucleonic 
matter this model gave bad results [TH [75] . 

2.1 Argonne NN interaction 

A typical NN potential contains a strong short-range repulsion, an intermediate- 
range attraction and a long-range part. 

The short-range part is usually treated phenomenologically, because for short 
distances and high momenta transferred the quarks inside nucleons start to play an 
important role which in principle can not be well included in an effective potential 
without using explicit subnuclear degrees of freedom. 

The intermediate-range attraction was initially phenomenologically described as 
a sum of Yukawa functions in the Reid|76j potential, or as the exchange of a scalar 
meson a as in the one boson exchange (OBE) potentials [77] . The long-range part is 
due to the one pion exchange (OPE) between nucleons. 

Other non-local terms like the linear and quadratic spin-orbit terms, as well as 
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relative angular momentum or other operators coming from two pion exchange or 
from nucleon A excitations, are usually included in modern interactions. 

The main contribution of NN potential at low energy can be derived by con- 
sidering that nucleons exchange a meson. The long-range part of Vij is know to 
be mediated by the pion vr, that is the lightest of all mesons. This is the common 
starting point of all recent and past NN interaction. The OPE potential is given by 

„,Tr JttNN ^ ^ /O 0^ 

= ^^ij^i ■ ' (2-2) 

with 

Xij = Y{m.„rij)ai ■ aj + T{m.„rij)Sij , (2.3) 

where 

^^^^ = 0.075 ± 0.002 (2.4) 
47r 

is the pion- nucleon coupling constant [78], m.,^ is the averaged pion mass, a and r are 

Pauli matrices acting on the spin or isospin of nucleons, Sij is the tensor operator 

Sij = 3{ai ■ fij){aj ■ rij) - ffi ■ Bj , (2.5) 

and the radial functions associated with the spin-spin and tensor parts are 

Yix) = — Cy(x), 



X 

T(x) = ( 1 + ^ + ^ ) y(x)er(x) . (2.6) 



The ^(x) are short-range cutoff functions defined by 

Cy(x)=CT(x) = l-e— (2.7) 

It is important to note that since T{m.,^r) ^ Y^niT^r) in the important region 
where r <2 fm, the OPE is dominated by the tensor part. 

The complete NN potential Vij is given by vfj + v(j, where v/^ contains all the 
other intermediate-range and short-range parts. All the parameters in the ^ short- 
range cutoff functions as well as other phenomenological vj^ parts are fitted on NN 
Nijmegen scattering data. 

The most recent interaction proposed by the Argonne group is called AV18[57J 
NN potential because it is written as a sum of 18 operators. This is a modern 
interaction fitting very well the scattering data of Nijmegen database. There are 
other simpler Argonne interactions called AVn'[58] that are simpler versions of AV18; 
they contain n < 18 operators, and the prime symbol indicates that such potentials 
are not just a simple truncation of AV18, but a reprojection, which preserves the 
isoscalar part of AV18 in all S and P partial waves as well as in the ^Di wave and 
its coupling to ^Si. 
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The Argonne NN potential between two nucleons i and j is written in the coor- 
dinate space as a sum of operators 

% = E M^^^)0^j ' (2-8) 

p=l,n 

where n is the number of operators depending of the potential, Vp{r) are radial 
functions, and rij is the inter-particle distance. 

The first eighth operators give the higher contribution to the NN interaction. 
The first six of them come from the long-range part of v^j OPE term. The last two 
terms depend on the velocity of nucleons and give the spin-orbit contribution. These 
eight operators are 

Oif^'^ = (1, ■ ffj, SijXij ■ Sij) X (1, Ti ■ fj) , (2.9) 
where Lij is the relative angular momentum of couple ij 

Lij = ^(^i - 5) X - ^i) ' (2-10) 

and Sij is the total spin of the pair 

Sij = \{oi + aj). (2.11) 

In the modern interactions these eight operators are the standard ones required to 
fit S and P wave data in both triplet and singlet isospin states. 

Following operators used in the Argonne AV14[67] and AV18 are 

O^f^^^ = (4,4(a, • a,), (4- • Sijf) X (l,r, • fj) . (2.12) 

These additional operators were introduced to better describe the phase shift in 
both P and D state |79j. The four L?' operators provide for differences between S 
and D waves, and P and F waves. In addition to S and L ■ S operators the (L • S)"^ 
provide an additional way of splitting states with different J values (for example 
■^-^^1,2,3 states). The effect of the quadratic angular momentum operators instead 
of and consequences were studied by Wiringa et al.[80]. However, the contribution 
of these operators is small respect to the total potential energy. 

The last four additional operators of the AV18 potential break charge indepen- 
dence and are given by 

QP=i5,i8 ^ j.^^.^ . S,,Tij,T! + rj , (2.13) 

where Tjj is the isotensor operator defined as 

= StIt^ - f ■ fj . (2.14) 

The Argonne AV18 potential has been used in a wide range of works about the 
ground-state [HI [351 (SHI [HD [821 IHS] and excited states [HI] of nuclei and nucleonic 
matter [85]. 
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2.2 Urbana and Illinois three-body forces 

The three-nucleon interaction (TNI) Vijk depends on the choice of the NN potential, 
but the final result with the total Hamiltonian should be independent of the choice. 
For this purpose several TNI form have been used, fitted in addition to the Argonne 
NN interactions, namely the Urbana IX|64j and the modern Illinois[63j forms. 

The Urbana and Illinois structure of the TNI interaction between three nucleons 
is written as a sum of several operators 

where each term represents the interaction given by two-pion exchange in the P and 
S wave, the three-pion exchange and a phenomenological spin-isospin independent 
term. The Urbana force only contains the first and the last terms, while the more 
sophisticated Illinois forces contain all of these operators. 

The first term was introduced by Fujita and Miyazawa[86] to describe the process 
where two pions are exchanged between nucleons with an intermediate excited A 
resonance as shown in Fig. 12. 1[ 









A 


7t 





Figure 2.1: Fujita-Miyazawa 27r-exchange in P wave term 



This diagram with three nucleons can be schematically described as 

(U) 



eye 



where the sum is cycled to include all possible permutations, and 

describes the pion exchange transition potential that scatters the |A^A^) to a \NA) 

state and viceversa. 

The OPE interaction has to be the same form of NN potential, that for the AV18 

is 

f 2 

~ 47r 3 ^ ' ^ ' 
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with X-j defined as Xij in Eq. 12.31 and the functions Y{x) and T{x) defined in Eq. 



12. 6[ The short-range cutoff functions are 

^y{r)=CT{r) = l-e-''\ (2.18) 



„2 



The Xij is the same interaction of pion-exchange transition potential entering in 
Eq. I2.16] that. with the inclusion of the appropriate spin-isospin operators, becomes 

VNN^NAi^j) = Vtr{ij)Xijai ■ S^Ti ■ , (2.19) 

and 

'"an^nnU^) = vtr{jk)XjkSj • (Tfcfj • ffc , (2.20) 

where a and r are the usual spin and isospin operators acting on nucleon, and 
ft are the spin-isospin transition operators between the to the A space, and S 
and T for the reverse transition. 
Using the Pauli identity 

■ at'' ■ B = -A- B -i-T- Ax B , (2.21) 
3 3' ^ ^ 

the Eq. 12.161 can be rewritten in the usual form 

^r2TT,PW _ .PWrP-J^'P^ — 

^ijk — ^27r '^ijk — 



eye 



1 

{Xij,Xjk} {n ■ Tj,Tj ■ Tk} + - [Xij,Xjk] [n ■ tj^tj ■ r^] 



(2.22) 



where the constant A2^ is fitted to reproduce the ground state of light nuclei and 
properties of nuclear matter, and is different for each model of three-body potential. 
In the Urbana forces the constant typically denoted by ^27r- 

The second term in Eq. 12.151 is due to the ttN scattering in the S wave as shown 
in Fig. [221 

The form of this operator is 

^Jfc"^^ = Mi^^O'fjf^ = ^ItT XI Z{mT,rij)Z{m^rjk)ai ■ hjC^k ■ rkjn ■ Tk , (2.23) 

eye 

with 

Z{x) = ^ [Y{x) - T{x)] . (2.24) 

This term is essentially a simplified form of the original Tucson model [87]. The 
structure is the same of anticommutator term in Eq. 12. 22) with a redefinition of radial 
functions and operators and the spin-isospin dependence reduces to a quadratic form. 
This term is not included in the Urbana forces. 

The three-pion exchange term consists in two diagrams where pions excite one 
or two nucleons into intermediate A-states as shown in Fig. 12.31 and [27 
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Figure 2.2: 27r-exchange in S wave term 







K 









Figure 2.3: 37r-exchange with a A-intermediate state 



The scattering potentials of the process of Fig. 12.31 where only one A is created 
can be written in the same way as described for the Fujita and Miyazawa term by 
considering the possible excitation because of the pion exchange: 



3tt,AR 



l,ijk 



E 

eye 



1 



[vNN^ANiij)vHjk)vlN-.NNiik) + j ^ k] , (2.25) 



where the symbol j ^ k means the exchange of nucleons j and k to consider all the 
possible permutations between nucleons. In the same way it is possible to describe 
the interaction of diagram in Fig. 12.41 that is 



3tt,AR 



1 



2,tjk — m/v) 

eye ^ ^ ^ ' ' 



2 l^NN~^NA 



iij)vAN^NAUk)vJfA-.NNiik)+j ^ k] , 

(2.26) 
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Figure 2.4: Svr-exchange with two A-intermediate states 



with the PauH identity given in Eq. 12.211 and the following 

ft.f = 2 

f t X f = -i-T 
3 

t-At-B = A-B + ia-AxB, (2.27) 

the two diagrams can be rewritten in terms of other operators as reported in Ref. 
|63| in a very complex way. 

The last phenomenological part of TNI is due to all neglected physical effects. 
This part has the role to compensate the overbinding and the large equilibrium 
density of nuclear matter given by other TNI operators. 

The is a spin-isospin independent operator given by 

V^, = AnOf^^ = ^rY. THm^n,)T\m^r,k) , (2.28) 

eye 

where T[x) is the same function defined in Eq. 12. 6i 

The A2tt term of Urbana UIX was originally fitted, in addition to the AV18, 
to reproduce the triton and alpha particle binding energy, while the Uq strength 
was adjusted to obtain the empirical equilibrium density of nuclear matter [M]. It 
has to be noted that the ground state of light nuclei can be exactly solved with 
few-body techniques, but the determination of the equation of state of symmetric 
nuclear matter can be evaluated only using many-body techniques that contain 
uncontrolled approximations. Although the AV18+UIX Hamiltonian gives very good 
results for the triton and alpha particle ground-state, a calculation of p-shell nuclei 
using GFMC revealed that this Hamiltonian cannot reproduce a binding energies in 
agreement with experimental data|83j. 

The more modern Illinois forces were fitted using GFMC to correctly reproduce 
all the 17 bound stated of 3<A<8 nuclei[63], and then they have been proved to 
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well describe the ground state of nuclei up to ^^C|35j. However, the Illinois TNI 
presents some problems if used to evaluate the equation of state of pure neutron 
matter, as it will be discussed in the following chapters, possibly because they were 
fitted without any information coming from nuclear matter calculations. 



Chapter 3 



Method 



The central problem in all many-body theories is the solution of the Schrodinger 
equation to evaluate the ground-state properties of a generic system. An exact ana- 
lytic solution is impossible to be found for most many-body Hamiltonians, and also 
the numerical brute- force solution is impossible in most interesting cases. 

The Diffusion Monte Carlo (DMC) method, or the Green's Function Monte Carlo 
(GFMC) method, that is very similar, are ways of exactly solving the many-body 
Schrodinger equation by means of a stochastic procedure. Although some attempts 
to develop such Quantum Monte Carlo (QMC) methods were just explored, the 
first application of a GFMC code is due to Kalos[SB], and in 1962 the ground state 
of three- and four-body nuclei using some simple central interactions was solved. 
The first DMC algorithm was implemented by Anderson in 1975 [89] to study some 
molecular ground-state, and the first simulation of uniform matter was carried out by 
Ceperley and Alder [90]. The application of the GFMC algorithm to nuclear Hamil- 
tonians containing spin and isospin operators was proposed in 1987 by Carlson [33]. 

The extension of the DMC to deal with nuclear Hamiltonians using a large 
number of nucleons was invented by Schmidt and Fantoni 18] with the Auxiliary 
Field Diffusion Monte Carlo (AFDMC) method, used to solve the ground state for 
a generic system given a nuclear Hamiltonian, where the interaction contains spin 
and isospin operators as described in chapter [2l 

In this chapter it will be review the DMC method to calculate the ground-state 
energy of a quantum system by starting from its Hamiltonian containing a kinetic 
term and a scalar interaction, then it will be explained the extension to the AFDMC. 



3.1 Diffusion Monte Carlo 

The DMC method [9T] 192] . projects out the ground state from a trial wave function, 
provided that it is not orthogonal to the true ground state. 
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The Schrodinger equation in imaginary time r = it is given by 

- ^i^{R,T) = Hi,{R,T) , (3.1) 

with R = {ri...rAr} is a configuration of the system with all its degrees of freedom. 
A formal solution written in integral form is 

'ijj{R,T)= j G{R,R',T)'il^{R',0)dR' , (3.2) 

where the kernel G{R, R' ,t) is the Green's function of the operator H + and can 
be expressed as the matrix element 

G{R,R',t) = (i?|e-(^-^»)n^') = j;e-(^"-^^)-<A:(i?)(/'„(i?') , (3.3) 

n 

where {4>n} is a complete set of eigenvectors of H. By considering a generic trial 
wave function ^jJT as a linear combination of {(pn}, it is possible to consider the 
evolution in imaginary time r, given by 

^P{R, r) = e-(^--^o)^V(-R, 0) , (3.4) 

where Eq is the exact ground state energy. In the limit r oo, tp{R,T) approaches 
to the lowest eigenstate "00 with the same symmetry of "0- The evolution can be done 
by solving then the integral of Eq. 13.21 in the limit of infinite imaginary time. 

Let is considered the non interacting Hamiltonian of N particles with mass m: 

H, = -Y^^Vl (3.5) 

i 

The Schrodinger equation 13.11 becomes a 3A^-dimensional diffusion equation, and 
by writing the Green's function of Eq. 13.31 in the momentum space by means of 
the Fourier transform, it is possible to show that Go is a Gaussian with variance 
proportional to r: 

m \ — 



This equation describes the Brownian diffusion of a set of particles with a dynamic 
governed by random collisions. This interpretation can be implemented by repre- 
senting ip hy a set of discrete sampling points that generally are called walkers. 
Therefore it can be defined as: 

i;{R) = Y,m-Rk)- (3.7) 
k 

The evolution for an imaginary time step At is then realized by 

ij{R,T + AT) = Y,GiR,Rk,AT). (3.8) 



3.2. IMPORTANCE SAMPLING 



19 



The result is a set of Gaussians that in the infinite hmit of imaginary time r represent 
a distribution of walkers according to the lowest state of Hamiltonian and can be 
used to calculate the ground state properties of the system. 

Now, let be considered a more realistic case with a generic Hamiltonian 

i 

By using the Trotter formula 

g(A+B)AT ^ ^-A^g-BAr^-A^ ^ o(At3) , (3.10) 

it is possible to break up the Green's function in two manageable parts, which can 
be written as: 

V(R) + V(r') . 

G{R, R', At) « e 2 ^^Go{R, R', At) . (3.11) 

This approximate expression is valid only in the limit At — > 0. Eg. 13.21 becomes: 

/ m \^ f m(R-R'f V{R)+V(r') a 

i,{R, t + At) = (^^^^) J e-^^e 2 ^^i^iR', T)dR' , (3.12) 

where the factor due to the interaction with the trial eigenvalue Et (that is a 
normalization factor) is the weight of the Green's function: 

The integral 13.121 can be solved in a Monte Carlo way by propagating the par- 
ticle coordinates by sampling a path according the diffusion term in the integral 
(Gaussians). The sampling of the weight term w is often realized with the branch- 
ing technique in which w gives the probability of a configuration to multiply at 
the next step according to the normalization. Computationally, this is implemented 
by weighting estimators according to and generating from each single walker a 
number of walkers 

n=[w + i\, (3.14) 

where ^ G [0; 1] is a random number and [x] means integral part of x. 

This is the most used technique to sample the weight term, but other techniques 
to sample w are possible but not interesting in this work. 



3.2 Importance Sampling 

The technique described in the above section is rather inefficient because the weight 
term of Eq. 13.131 could suffer of very large fluctuations, for example when a particle 
moves close to another one. 
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A big improvement to overcome this problem is possible by means of importance 
sampling techniques. The Green's function used for propagation can be modified as 
following: 

G{R, R', At) = ^ytUg{R, R' , Ar) . (3.15) 

The so called importance function ipj in the above equation is often (but not 
necessarily) the same of that used for the projection. In this case the walker's dis- 
tribution is sampled according the function 

f{R,r) = i^i{R)ij{R,T). (3.16) 

By inserting Eq. 13.151 in 13.121 and expanding near R, the integration gives an 
additional term in Gq and in the weight to keep the normalization of G. The diffusion 
is now guided by an additional drift term 

Go{R, R', At) ^ Go{R, R' , Ar) = (^^^^) e , (3.17) 



and the weight 13.131 becomes 



El(R)+El(R') 



-Et At 



tt) — > t() = e V / J (3.18) 

where 

is the local energy of the system. 

It is possible to show that if one multiply the Eq. l3.1l bv the importance function 
ipl, the Schrodinger equation tunes into a Fokker-Plank equation: 



^fiR,r) = -^V'fiR,r) + ^V 
OT Zm 2m 



+ EL{R)fiR,T), (3.20) 



and it can be verified that the corresponding Green's function is given by the two 
terms of Eq. [3171 and [3T51 



3.3 Auxiliary Field Diffusion Monte Carlo 

In the case of nuclear Hamiltonians the potential V contains quadratic spin and 
isospin and tensorial operators, so the many body wave function cannot be written 
as a product of single particle spin-isospin states. 

For example, let consider the generic quadratic spin operator ffi ■ aj where the 
a are Pauli's matrices operating on particles. It is possible to write 



fj,: • (To 



0-,; 



CT,- 



2Pj 



1, 



(3.21) 
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where P-j interchanges two spins, and this means that the wave function of each 
spin-pair must contain both components in the triplet and singlet spin-state|93t [94] . 
By considering all possible nucleon pairs in the systems, the number of possible 
spin-states grows exponentially with the number of nucleons. Another way to un- 
derstand this problem is simply connected to the fact that such operators cannot 
be diagonalized in a single particle spin-base. These problems are obviously related 
also to quadratic isospin-operators. 

In order to perform a DMC calculation with standard nuclear Hamiltonians, it 
is then necessary to sum over all possible single particle spin-isospin states of the 
system to build the trial wave function used for propagation. This is the standard 
approach in GFMC calculations for nuclear systems. 

The idea of AFDMC is to rewrite the Green's function in order to reduce the 
quadratic dependence by spin and isospin operators to linear by using the Hubbard- 
Stratonovich transformation. 

Let us consider the first six operators of the NN potential of the form of Eq. 12.91 
The operators can be recast in a more convenient form 

6 

^ = E E Mr^j)0^^''Hhj) = Vsi + VsD = 

(3.22) 

where Latin indices label nucleons, Greek indices stay for Cartesian components, 
and 

Vsi = Y.vi{nj) (3.23) 

is the spin-isospin independent part of the interaction. The 3A by 3A matrices A^'^^ 
and A^'^'^^ , and the A by A matrix A'^'^^ contain the interaction between nucleons of 
other terms: 

= '^4(r-ij)(5c,/3 + ^^6(rii)(3fgfg -Sap), 

At]=V2{n,). (3.24) 

The matrices A are zero along the diagonal (when i = j), in order to avoid self in- 
teraction, and are real and symmetric, having then real eigenvalues and eigenvectors 
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given by: 

E^S^S = ^^VS- (3-25) 
j 

The matrices A*^*^) and A^'^'^^ have n = 1...3A eigenvalues and eigenvectors, while 
A^'^^ has n = 1...A. One can then define a new set of operators written in terms of 
eigenvectors of matrices A: 

Ol^2 = Y.^^ai^n^■ (3-26) 

The potential 13.221 becomes 

3A 3 3A 3 A 

n=l Q=l n=l Q=l n=l 

Now it is possible to reduce the spin-isospin dependence of the operators from 
quadratic to linear by means of the Hubbard-Stratonovich transformation. Given a 
generic operator O and a parameter A: 



e-k^o^ = J- I dxe-^+v^-O _ (3 28) 



.2 



With the Hubbard-Stratonovich transformation it is possible to linearize the quadratic 
form of the new set of operators On defined in I3.26[ Then the Green's function be- 
comes: 

X TT / dx„e-^e^-^"^"^"^" , (3.29) 

n=i V 2vr J 

the On contain all the 3A On"^ operators, the S)A oi^J^ operators, and the 3^ 

(t) 

Ona operators. The newly introduced variables x„, named auxiliary fields, have 
to be sampled in order to evaluate the integral in Eq. 13.291 The linear form of 
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this operators gives the possibihty to write the trial wave function as a product of 
single particle spin-isospin functions. In fact the effect of the new operators On in 
the Green's function consists in a rotation of the spin-isospin of nucleons without 
mixing the spin-isospin state of nucleon pairs. 

The sampling of auxiliary fields to perform the integral in Eq. 13.281 eventually 
gives the same effect as the propagator with quadratic spin-isospin operators acting 
on a trial wave function containing all the possible good spin-isospin states. The 
effect of the Hubbard-Stratonovich is then to reduce the number of terms in the 
trial wave function from exponential to linear. The price to pay is an additional 
computational cost due to the diagonalization of A matrices and the sampling of 
the integral over auxiliary fields. 

Sampling of auxiliary fields can be achieved in several ways. The more intuitive 
one, in the spirit of Monte Carlo sampling is to consider the Gaussian in the integral 
13.281 as a probability distribution. The sampled values are then used to determine 
the action of the operators on the spin-isospin part of the wave function. This is 
done exactly as in the diffusion process by noting that a each field value doesn't 
depend by the old value. Other techniques to solve the integral are possible (i.e. 
with the three-point Gaussian quadrature [95|) but they will not be considered in 
this work. 

The introduced auxiliary fields might also be given a physical interpretation. 
In fact, they might be seen as a sort of meson field that generate the interaction 
between nucleons. By considering the simple one-pion exchange picture between two 
nucleons, the operators in the Green's function 13.291 represent exactly the interaction 
between a pion field with a nucleon. In fact the pion couple to nucleon with a term 
of the form: 



where tTq, are the vr and vr fields. 

By considering fixed nucleon coordinates, the faster pion degrees of freedom 
are integrated out giving the NN interaction as usual. In this way the Hubbard- 
Stratonovich is a reverse way of re-introducing the pion fields responsible of NN 
interaction. The fact that in sampling of integral 13.281 each value of auxiliary fields 
does not depend on the previous one means that a pion fields does not have any 
correlation during the diffusion of slower nucleons, and the meson's kinetic energy 
can be included in counter terms of the interaction. This is exactly what is com- 
monly done by integrating out the meson field to obtain a meson independent NN 
interaction. 

The importance sampling in the Hubbard-Stratonovich transformation that ro- 
tate nucleon spinors can also be included. For auxiliary fields importance sampling 
is achieved by "guiding" the rotation given by each 0„ operator. More precisely one 
can consider the following identity: 



(3.30) 



x: 



,2 



x: 



.2 



f + V^XniOn) + V-^uAt {On " (0„)) , (3.31) 
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where the mixed expectation value of the operator is calculated in the old spin- 
isospin configuration: 

This can be implemented by shifting the Gaussian used to sample auxiliary fields, 
and considering the extra terms in the weight for branching 

where 



X = yJ-\n/^T{On) ■ (3.34) 

The additional weight term in Eg. l3.33l to be added to Eg. 13.291 can also be included 
as a local potential, so it becomes 

e ('^/i^-s) . (3.35) 

By combining the diffusion, the rotation and all the additional factors it is possible 
to obtain two equivalent algorithm. The first explicitly contains the drift correction 
included in the Gaussian used for the diffusion process and the extra terms given 
by the guidance of rotations: 



where the drift term d can be arbitrarily fixed (for example, one could choose to 
guide the diffusion according to the modulus or with the real part of the trial wave 
function), by only correcting the factor ^^^jj^g^^ ■ 

The second is calculated in a simpler 'local energy' scheme, as described in 
section 13.21 and derived in Eq. I3.18t 



exp 



v'm\RS)\ mv\RS) \ . 

2m \{MRS)\ mRS) ) 



{MR'smi\RS)\ 

{^Iji\RS)\{^i\R'S)\ 



(3.37) 



where in the last expression the choice of the drift term is taken to be 



(3,38) 

but the generalization to a general drift term is trivial. 

If the drift term is the same, the two algorithms are equivalent and should sample 
the same Green's function. 
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3.4 Spin-orbit propagator for neutrons 

As described in Sec. 12. li the nuclear Hamiltonian contains also a spin-orbit term that 
has to be treated separately from other terms described in the previous section. 
The spin-orbit potential is defined 



VLs{hj) = VLs{rij)Lij ■ Sij = VLs{rij){L ■ S)ij 



(3.39) 



Because the spin-orbit is a non-local operator the corresponding Green's function is 
not trivial to be derived. 

One way to evaluate a L- S propagator is to consider the first order of expansion 
in At 



e 2 



VLs(rij)iL-S)ijA7 



1 - ^VLs{rij){L ■ S)ijAT 



(3.40) 



and acting with this on the free propagator Go[93]. The derivative terms of the 
above expression give 



-» -* 711 

(Vi - Vj)GoiR,R') = -j^i^n - Af,)GoiR,R') , (3-41) 



where Ar^ 
Then 



r,- — r-. 



{L ■ S),,Go{R,R') = -J^J^in - rj) X {An - Ar,) • {di + a,)Go(i?, i?') , (3.42) 

and after the exponentiation and the multiplication of VLs{f)j the following propa- 
gator is obtained: 



Pls = exp 



Y.J^''Ls{rij){L-S)ij 

i+3 



(3.43) 



This propagator must be corrected for counter terms arising at order At from the 
approximation implied in Eq. 13.401 Note that in Eq. 13.421 there are terms depending 
on V At in the integration for the evolution of wave function in Eq. 13. 2t 



r X r' = r X (r + Ar) = r x Ar ~ V At , 



(3.44) 



and one expects that other spurious terms linear in Ar are intrinsically contained 
in Pls- 

In order to see the effect of these additional terms the trial wave function is 
expanded to the first order: 



^{R') ~ *(i?) + ^Vi^'(i?)( 



r,; - ri) + 



(3.45) 
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with i? = "fi, ..,r/v Then the Pls is also expanded, giving 



Pi5 ~ 1 - ^ 



4i 



EE 

ij^j k^i 



m'^VLs{rij)vLs{rki) 



(Eij X fij) ■ Afi(Sfc/ X fki) ■ Afk + ... 



where 



ai + a 



(3.46) 
(3.47) 



and it was used the relation a ■ {h x c) = c ■ {a x h). 
By inserting Pls with the Go in the Eq. 13.21 gives 



^'(i?,T + Ar) = j dR'Go{R,R',A 

m'^VLs{rij)vLs{rkl) 



i-E 



mvLs(rij) ^ 



EE 

ij^j k^i 



32^4 



X Tij) • Afi{T,ki X rfe/) • Affc + 



^(i?,T) + ^Vi^(i?,r)( 



r,; - + ... 



(3.48) 



The integration in the variables R' shows the presence of extra terms. Note that 
terms linear in Ar and quadratic with different components (like AfiAfj with i ^ j) 
integrate to zero, while quadratic terms in the same components give ^^Ar/m. The 
integration of Eq. 13.481 of the part with quadratic Af coming from linear terms of 
both Pls and the wave function gives 



m '■ — ' 4z ft^ 



-ATY,VLsirij){l ■ S)ij^{R) , 



(3.49) 



that is the spin-orbit contribution of the propagator. 

The integration of the quadratic terms in Ar coming from Pls gives the addi- 
tional spurious terms to the propagator. These extra corrections have the structure 
of a two- and three-body terms: 



AT{y2 + V3) = ^"^^^ X] X] X] '^^^sk<'ij)VLsiTih){X^ij X Tij) ■ (Sjfe X fjfc) 

At— 2 X] H ^Lsirij)vLsirik)[rij ■ fik^ij ■ E^fc - S^^ • 



i jf^i k^i 



(3.50) 
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The term with j = k is an extra two-body contribution 



add 



E 

i<j 



[2 + ai- Uj - Ui ■ fijaj ■ fij] 



(3.51) 



while that with j ^ k gives an additional three-body correction 

mfij ■ rikVLs{rij)vLs{nk) 



■EE 

i<j<k eye 



[fij ■ fik{2 + Gi ■ Gj + ai- Gk 



-Gj ■ Gk) - Gi ■ TijGj ■ Tik - Gk ' VijGi ■ Tik - Gk ' TijGj ■ Tik] 



(3.52) 



Note that both the two- and the three-body additional terms contain a quadratic 
form of spin operators, and in order to be sampled they need some additional 
Hubbard-Stratonovich variables. 

An alternate way to include the additional counter terms given by this spin-orbit 
propagator is to consider a different propagator of the form 



exp 



X exp 



El m A ^ 1 ^ 



(3.53) 



This different propagator contains the quadratic form of spin operators and need 
additional Hubbard-Stratonovich fields to be applied. However, by expanding the 
two exponential to keep terms linear in Ar, the above expression gives 



^ + 2^ Ji^'"Ls{rij)[rij X Arij] • Gi 



1 m? 



-^^VLs{rij)vLs{rklWij X Arij) ■ Gi{fkl x Arkl) ■ ^k + 

2 



X 



+ ^ E E J^-^'"Ls{rij)vLs{rki){rij x Arjj) • Gi(fki x Aru) ■ ^k 



(3.54) 



the quadratic forms of Ar giving the extra terms after the integration cancel each 
selves and others are of higher order in Ar. The two forms of spin-orbit propagator 
are equivalent to the first order in Ar. 

An important remark concerns the spin-orbit propagator in presence of isospin- 
operators. For neutrons the r^ • tj operators is a constant, but the inclusion of these 
operators to the L-S is hard. The additional counter terms without isospin-operators 
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have a three-body structure, but the spin-operators appearing there are quadratic, 
and the Hubbard-Stratonovich transformation can be used safely. In presence of 
isospin-operators, the counter terms contain cubic spin-isospin operators. This fact 
prevents the straightforward use of the Hubbard-Stratonovich transformation. 

3.5 Three-body propagator for neutrons 

For pure neutron systems the three-body forces in the Urbana or Ihinois form can be 
rewritten as a two-body term, and then safely included in the Hubbard-Stratonovich 
transformation as described in the above section. 

Urbana and Illinois three-nucleon potentials can be written in the form 



T/ _ APWr,'^TT,PW ^SWr,'^'rr,SW .ARr,^TT,AR j_ A 



Let us examine the structure of each of these operators. As explained in Sec. 12.21 
the one pion exchange interaction at the base of the three-body force is the same 
of that used in the NN interaction, implying that the OPE operator X°? defined in 
Eg. 12.31 has the same algebric structure as the spin dependent part of the two-body 
potential Vsd oi Eq. I3.22[ Therefore, it can be expressed in a similar way, in terms 
of a 3A^ by 3A^ matrix Xi^-jj^ as 



with Xia-jp = if i = j- Notice that this matrix is symmetric under Cartesian 
component interchange a ^ [3 and under particle label interchange i ^ j and is 
also fully symmetric Xj^-i^ = Xia-j/s- 

For neutrons the commutator term in Ofj^'^^ vanish, and the anticommutator 
can be reduced in a quadratic form of the spin operators. Thus the operator of Eq. 
12.221 reduces to 




(3.56) 




(3.57) 



i<j<k 



i<j 



where 




(3.58) 



k 



ihe iucson b wave component C^^.^ ( 
operators and it is convenient to write it as: 



of the Eq. 12.231 is quadratic in the spin- 




(3.59) 



i<j<k 



i<j 



k 



where 



Z{x) 



^ [Y{x) - T{x)] 




(3.60) 
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The 3-pion exchange O^J^^^ terms have a central part and a spin dependent 
part. 



i<j<k 



The central part is given by 



400 2 y 



which is 200/9 times the trace of X^. The spin dependent part is 



(3.61) 



(3.62) 



200 



(3.63) 



and one can just write out the 4 nonzero terms for each combination of ai and ffj. 
Finally the spin independent terms can be written using just pair sums as 



E ^^jk — + - y^(G 



B.\2 
i J 1 



(3.64) 



with 



i<j<k 



y^T^(?T7-^rifc) 



^3) ■ 



(3.65) 



i<j 



Then the spin dependent part of the three-body interaction can be easily in- 
cluded in the matrix Ai^-jp by 



jk 



200 



I aARy2 y c c 

and the central term to be included in the propagator is given by 



(3.66) 



VsiiR) ^ Vsi{R) + Ar\gI^ + -Y,iG. 



R\2 



,400 



"T ^Stt 2-^ ^ia;jf3^ia;jl3 ■ 



i<j 



(3.67) 



The coefficients of three-body forces Urbana and Illinois can be found in Ref. 



ESI 
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3.6 The constrained-path and the fixed-phase approxi- 
mation 

As described in the above sections, the DMC project out the ground-state of a 
given Hamiltonian as a distribution of points in the configuration space. On the 
other hands for the diffusion interpretation to be vahd, the trial wave function must 
always be positive definite, since it represents a population density [96j. Thus the use 
of DMC is restricted to a class of problem where the trial wave function is always 
positive such as for a many-Boson system in the ground-state. 

For Fermionic system the DMC algorithm can be used by artificially splitting 
the space in regions corresponding to positive and negative regions of the trial wave 
function. It is possible to define a nodal surface where the trial wave function is zero 
and during the diffusion process a walker acrossing the nodal surface is dropped; this 
approximate algorithm is called fixed- node |97| 198]. It can be proved that it always 
provides an upper bound to the true Fermionic ground-state. 

Other difficult techniques to overcome the Fermion sign problem exist, like the 
transient estimators analysis [90] or other alternative methods to control the sign- 
problem [99l[T00l[l0T]. However, at present, none of these methods is efficient enough 
to compute physical properties with a sufficient accuracy (with the noticeable ex- 
ception of the electron gas). 

In the case of nuclear Hamiltonians, or for each problem where the trial wave 
function must be complex, the constrained-path approximation [T02| 11031 I104j is 
usually apply to avoid the Fermion sign problem. The constrained-path method was 
originally proposed by Zhang et al. as a generalization of the fixed-node approxima- 
tion to complex wave functions. 

If the overlap between walkers and the trial-wave function is complex (as in our 
case), the usual sign problem becomes a phase problem. It is possible to deal with it 
by constraining the path of walkers to regions where the real part of the overlap with 
the trial wave function is positive. The first AFDMC calculations were performed 
by applying this constrain [18j. 

Let be considered a complex importance function. In order to keep real the space 
coordinates of the system, the drift term of Eq. 13.171 that gives the importance 
sampling has to be a real term. In the case of the constrained-path approximation, 
a good choice for the drift is ^ 



and to eliminate the decay of the signal-to-noise ratio it is possible to impose the 
constrained-path approximation, that is realized by requiring that the real part of 
the overlap of each walker with the trial wave function must keep the same sign. 

^It should be possible to consider a different drift, but this form was found to give better results 
and was employed in the past AFDMC calculations 



d 




(3.68) 
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Thus, one has to impose that 

Re[iPi{R')] 
Re[iPj{R)] 



>0, (3.69) 



where R and R' denote the coordinates of the system after and before the diffusion 
of a time-step. If this condition is violated, the walker is dropped. 

An alternate way to control the sign problem is the fixed-phase approximation, 
that was proposed for systems whose Hamiltonian contains a magnetic field [XQ5j 
following earlier explorations by Carlson for nuclear systems [33]. 

Let us start with the same assumption that spacial coordinates of the system 
must be real, and let then consider the following drift term: 

Yltim (3 70) 

iV'/(^)i ■ ^^-^"^ 

With this choice the weight for branching becomes 

v^\MR)\ , VtPiiR) 



exp 



+ ^-^ Ar 



2m \MR)\ ^l{R) 



mR)\ MR') (^^,. 

\i;i{R')\ MR) ■ ^ ' 



Note that in the above expression there is the usual importance sampling term as 
in Eq. 13.151 and an additional term that corrects the choice of the drift term. 
A generic complex wave function can be written as 

^(i2) = |V(i?)|e^'^(«), (3.72) 

where (j){R) is the phase of ip{R), then the term appearing in Eq. 13.711 can be 
rewritten as 

\MR')\MR) ■ ^ ^ 

The fixed-phase approximation constrains the walkers to have the same phase 
as the importance function tpj. It can be applied by keeping the real part of the last 
term. In the same way as adopted for the importance sampling, after the expansion 
and the integration of the Green's Function to keep its normalization fixed, one has 
an additional term in the Green's function due to the phase, that must be included 
in the weight 

l-£(V<A)'Ar, (3.74) 

which can efficiently be included in the weight by keeping the real part of the kinetic 
energy. In fact: 

p WKfi) "^'IMR)] .v7^.m^2 

Re , = I . (V0CR . (3.75) 

ipiiR) \tpi[R)\ 

Then the real part of the kinetic energy includes the additional weight term given 
by the fixed-phase approximation. 
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Both the constrained-path and the fixed-phase are approximations to deal with 
the Fermion sign-problem and in principle should be equivalent if the importance 
function is close to the real ground-state of the system. 

It is important to note that Carlson et at. showed that with the constrained-path 
approximation the DMC algorithm does not necessarily give an upper bound in the 
calculation of energy [106]. This was also observed by Wiringa et al. in some nuclear 
simulations using the GFMC technique [82 j. 

Unfortunately it cannot be proved that the fixed-phase approximation gives or 
not an upper bound to the real energy. However, it was observed by Francesco 
Pederiva |107j that in some particular case the fixed-phase gives systematically en- 
ergies higher than the fixed- node energies (that always is an upper bound). The 
system studied was a quantum dot in an open-shell configuration; in that case the 
trial wave function is generally complex but can be also written as a real function 
with properly linear combination of the single-particle orbitals (for more information 
on a generic trial wave function, see the following section). 

3.7 Trial wave function 

The trial wave function used as the importance and projection function for the 
AFDMC algorithm has the following form: 

^t{R, S) = <^s{R)^a{R, S) , (3.76) 

where R = (ri, . . . , rV) are the Cartesian coordinates and S = {si, . . . , sn) are the 
spin and isospin coordinates of the system. The spin-isospin assignments Sj consist 
in giving the four-spinor components, namely 

/ ai\ 

_ bi 

Si = 

Ci 

\d^ J 

where Oj, bi, Ci and di are complex numbers, and the {\p t),p i), \n ]), \n |)} is the 
proton- up, proton-down, neutron-up and neutron-down base. 

The Jastrow correlation function is symmetric under the exchange of two parti- 
cles. Its role is to include the inter-particle correlation to the trial wave function for 
short distances. The generic form for the Jastrow is 

'^siR) = l{fjirij), (3.77) 

i<j 

where the function fj has been taken as the scalar component of the Fermi Hyper 
Netted Chain in the Single Operator Chain approximation (FHNC/SOC) correla- 
tion operator Fij which minimizes the energy per particle of nuclear or neutron 
matter ^lOHt I109j at the correct density depending to the system to be studied. 



= ai\p T) + bi\p I) + Ci\n T) + di\n j) , 
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The Jastrow part of the trial wave function in the AFDMC case only has the role 
of reducing the overlap of nucleons, therefore reducing the energy variance. Since it 
does not change the phase of the wave function, it does not influence the computed 
energy value in projection methods. 

The antisymmetric part of the trial wave function depends on the system to 
be studied; this function is generally given by the ground-state of non-interacting 
Fermions, that is usually written as a Slater determinant 



<^a{R,S)=A 



N 



\f ii 



Det{(pairi, Si)} , 



(3.78) 



where a is the set of quantum numbers of single-particle orbitals (pa depending on 

the system to be studied. 

For both nuclei and neutron drops orbitals are labelled with the set of quantum 
numbers a = {n,j,mj}: 



<t>a{ri, Si) = Rnjin) [Yi^jnii^)Xs 



(3.79) 



where i?„ j is a radial function (whose determination will be described in each 
specific case), i1,m( are spherical harmonics and Xs,ma are spinors in the usually 
proton-neutron-up-down base. The angular functions are coupled to spinors using 
the Clebsh-Gordan to have orbitals in the {n,j,mj} base according to the usual 
shell-model classification of the nuclear single-particle spectrum. 

In the case of nuclei or neutron drops, a summation of different determinants 
might be needed to build a trial wave function with the same symmetry of the 
ground-state of the nucleus considered. Because the AFDMC projects out the lower 
energy state not orthogonal to the starting trial wave function, it is possible to study 
a state with given symmetry imposing to the trial wave function the total angular 
momentum J experimentally observed. This can be achieved by taking a summation 
over a different set of determinants, thus 



D{R, S) 



(3.80) 



where the c/j coefficients are determined in order to have the eigenstate of total 
angular momentum J = ji + ... + jn ■ 

For a nuclear or neutron matter calculation the antisymmetric part is the ground 
state of the Fermi gas, built over a set of plane waves. The infinite uniform system 
is simulated with N nucleons in a cubic periodic box of volume L^. The momentum 
vectors in this box are 



27r 

T 



nn 



(3.81) 



where a labels the quantum state and Ux-, Uy and are integer numbers describing 
the state. The single-particle orbitals are given by 

0a(r;,Si) = e*^°-^^»Xs,m.(Sj)- 



(3.82) 
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The system has a sheh structure that must be closed, in order to meet the 
requirement of homogeneity and isotropy, then the total number of Fermions in a 
particular spin-isospin configuration must be 1, 7, 19, 27, 33,... It is also possible to 
modify the trial wave function to deal with an arbitrary number of Fermions and 
keeping zero the total momentum of the system [TTU] but this was not applied to the 
calculations discussed in this work. 

A particular phase of the system can also be simulated, for example to study the 
gap of a superfluid system. This is possible by including in the Slater determinant 
some pairing functions as done in several works about superfluid Fermi gases [TTT| 
11121 1113j . or by pairing Fermions with a Pfaffian[50 l I114j . 

In neutron matter the paired orbitals used in the Pfaffian are taken to be 

Mru , Si, Sj) = ^ ^e*^^-^^-C(5i, sj) (3.83) 

where the (,{si,Sj) is a spin function of two neutrons that correctly describe the 
spin-coupling, that in this case is in the ^S'o-channel, and vr and are BCS 
coefficients to be determined variationally. 



Chapter 4 



Results: nuclei and neutron drops 



In this and in the following chapters I will give a description of the results obtained 
by applying the AFDMC method to the study of nuclear and neutron matter and 
nuclei. 

This chapter is devoted to the analysis of finite systems, nuclei and neutron 
drops. The main concern is the comparison of AFDMC results with those computed 
by other methods. This first step is necessary to provide a benchmark test to assure 
the high accuracy of AFDMC calculations in dealing with Hamiltonians containing 
tensor forces. Some results on heavier nuclei, which can be tackled only by AFDMC, 
will also be given. 

The accuracy of AFDMC in calculating properties of neutron drops with a more 
sophisticated Hamiltonian with a spin-orbit and a TNI interaction in addition to 
the tensor force, and a comparison with the GFMC will be shown. The role of the 
TNI interaction for heavier neutron drop with an increasing of the density will be 
discussed. 

4.1 Open and closed shell nuclei 

As shown in chapter [2] the nuclear Hamiltonians are more complicated than those 
commonly used in condensed-matter problems. The increasing sophistication of NN 
and TNI interactions was not followed by a parallel development of a computational 
technique providing enough accuracy to predict properties of medium heavy nuclear 
systems, and progresses in accurate ground-state calculations has been quite slow. 

GFMC calculations was originally applied to the ground-state calculation of the 
alpha particle in 1987. The most recent progress was achieved in 2005 with the 
application to ^^C. In table 14. ll a chronological sequence of GFMC computations for 
heavier and heavier nuclei is displayed. 

Unfortunately the GFMC, that at present is probably the most accurate method 
applicable up to A=12 nuclei, will not be useful for nuclei heavier than A>12 at least 
for the future years because the huge power computing resources needed ^llGl 1117] . 
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nucleus 


year 


Ref. 


^He 


1987 


IMIiiSj 


^He 


1993 


[93] 


^Li and ^He 


1995 


m 


^Li and ^He 


1997 


[83] 


A=8 nuclei 


2000 


[82] 


A=9,10 nuclei 


2002 


m] 


12c 


2005 


[3S] 



Table 4.1: The ground-state calculations performed with the GFMC technique of 
nuclei up to the present. All the reported calculations were performed with some 
state-dependent nuclear Hamiltonian, and in the second column the year of publi- 
cation of corresponding results is reported. 



At present, the major goal in nuclear physics is the development of a many-body 
technique as accurate as the GFMC but applicable to medium nuclei. 

In this section it will be shown the capability of the AFDMC technique to deal 
with heavier nuclei problems by keeping the same accuracy of few-body techniques. 

4.1.1 Comparison between CP-AFDMC, FP-AFDMC and other 
techniques 

The possibility of applying the AFDMC with the constrained-path approximation 
(CP-AFDMC) to nuclei was originally explored by Schmidt et al. Early calculations 
yielded good results for nuclear Hamiltonians with a f4-like NN interaction [118], but 
in the presence of tensorial forces there was a very poor agreement with estimates 
obtained with other methods, as can be seen in table 14.21 where the CP-AFDMC 
results are taken from |119j . 



nucleus 


potential 


Eafdmc 


leather 


^He 


UV14 


-20.7(3) 




^He 


AV8' 


-17.9(1) 


-22.8(2) 


160 


UV14 


-84(2) 


-82.40 


160 


AV14 


-57(3) 


-84.0 



Table 4.2: The CP-AFDMC ground-state energy in MeV for the indicated potentials 
truncated to the vq level compared with other results [TT9] . The ^He energy was 
calculated with the GFMC using the Argonne AV8' [T20] . while the binding of the 
oxygen were performed with the CVMC[45] and with the CBF |12l| I122j using the 
Urbana t>iA|79j and the Argonne AV14[_67] interactions. 
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Both the alpha particle and the oxygen ground-state were calculated using a 
different NN interactions, namely the Urbana UV14[79] and the Argonne AV14[67] 
truncated to include only the first six operators. The alpha particle ground-state 
was solved with the GFMC method [120j . that cannot be used for oxygen. For the 
latter the only microscopic calculations available are the Cluster Variational Monte 
Carlo (CVMC)[15] and the Correlated Basis Functions fCBFl fml [T22] . 

The effect of the constrained-path approximation is apparently very different for 
the two interactions: for the AV14 interaction the CP-AFDMC does not reproduce 
well the binding energy of the alpha particle with a difference of about 1.2 MeV per 
nucleon higher than the GFMC estimate. In the case of oxygen with the UV14 the 
CP-AFDMC energy is roughly on par with the CBF calculation, while the result 
with the AV14 is very far from that of other variational calculations. It has to 
be noted that in the AV14 interaction the strength of the isoscalar-tensor and the 
isovector-tensor operators is higher than in UV14. Hence, the CP-AFDMC seems 
to well project the ground-state of a t;4-like interaction but not the isovector-tensor 
component of the potential. This fact suggests that such component induce a strong 
change in the nodal/phase structure of the many-nucleons wave function. 

These calculations were then repeated using the fixed-phase rather then the 
constrained-path approximation. The trial wave functions were modified with re- 
spect to the calculations of Ref. [119] . Radial orbitals were numerically computed in 
a self-consistent Hartree-Fock potential well generated using Skyrme forces [30 1 1123j 
with parameters fitted for light nuclei [12^ . The fixed-phase (FP) approximation 
turns out to be quite more effective than the constrained-path. 

The FP-AFDMC energy for alpha particle compares very well with the estimates 
from other methods, within 1% as can be seen in table HT3l 



method 


E 


FP-AFDMC 

GFMC 

EIHH 


-27.13(10) 
-26.93(1) 
-26.85(2) 



Table 4.3: The FP-AFDMC ground-state energy of the alpha particle in MeV for the 
AV6' interaction compared with other results [U]. The GFMC result is that reported 
in Ref. [58] after the subtraction of the Coulomb contribution of 0.7 MeV |125] . while 
the EIHH result is for the pure AV6' interaction [126j. 

The AV6' interaction employed [127] is a simplification of the full AV18 as ex- 
plained in section [2T1 FP-AFDMC energy is compared with the GFMC result of Ref. 
[58| after the subtraction of the Coulomb term that is of about 0.7MeV |125j . The 
other result were obtained with the Effective Interaction Hyper spherical Harmonic 
(EIHH)[T28l [129] that does not contain the Coulomb interaction [l26]. Both the 
GFMC and EIHH methods were proved to be in an excellent agreement with other 
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few-body techniques [25]. and the fact that FP-AFDMC is in par with this methods 
is the first important test to show its accuracy for hght nuclei calculation [47j. 

In order to have another comparison with GFMC, FP-AFDMC was used to 
calculate the ground-state energy of ®He[17] (note that the A=8 nucleus is just 
beyond the capability of few-body techniques). The AV6' interaction was employed. 
The ^He is an open-shell nucleus and in order to have a trial wave function which 
is an eigenstate of total angular momentum J = a linear combination of different 
Slater determinants is needed. Results are reported in table 14.41 



method 


E 


FP-AFDMC 
GFMC 


-23.6(5) 
-23.6(1) 



Table 4.4: The calculated FP-AFDMC ground-state energy in MeV for **He 
nucleus [U] with AV6' interaction compared with the GFMC result, whose result 
was taken from Ref. [58] and the Coulomb term of 0.7 MeV |125] were subtracted. 

It should be noted that ^He is an open shell nucleus. This fact does not seem to 
affect the agreement between FP-AFDMC and the GFMC, which is excellent. One 
more the calculation was performed using the AV6' interaction, after the subtraction 
of the Coulomb term from the GFMC calculation. An important further step has 
been the computation of the ground-state energy of ^^O. This calculation is way 
beyond the limits that presently affect GFMC, and comparison is possible only with 
variational calculations (either Monte Carlo or CBF based). FP-AFDMC overcomes 
the problems observed in previous CP-AFDMC calculations. 

The interaction used is the AV14 cut to the ue-level. Results are reported in table 
WM As it can be seen, the FP-AFDMC estimates is -90.8(l)MeV against -57(3)MeV 
of the CP-AFDMC. The effect of the constrained-path in this case is apparently 
catastrophic. The FP-AFDMC result is about 7% lower than other variational es- 
timates, consistent with the average difference between VMC and GFMC in light 
nuclei [Hi]. 



method 


E 


FP-AFDMC 

VCMC 

CBF 


-90.8(1) 
-83.2 
-84.0 



Table 4.5: The calculated FP-AFDMC ground-state energy in MeV for ^^O 
nucleus [17] with AV14 interaction cut to the ve-level. The VCMC [35] and CBF [122] 
results of the corresponding interaction are also reported. 
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The results on O show that AFDMC within the fixed-phase approximation 
can definitely deal with realistic nuclear Hamiltonians also for medium nuclei. 

4.1.2 Application of FP-AFDMC up to ^"Ca 

The capability of FP-AFDMC to treat medium-heavy nuclei can be exploited to 
understand some features of the interactions commonly used in literature. One im- 
portant open question is the portability of potentials fitted to reproduce light nuclei 
for heavier systems. How important the detailed operatorial structure does become? 
Can reprojections of sophisticated AV18 potential like AV6' be accurate enough to 
reproduce the main physical features of medium heavy nuclei? To this end it could 
be interesting to analyze the contribution to the total energy of nuclei with respect 
to the experimental data. 

The ground-state energy estimates of ^^O and ^^Ca, in addition to those of alpha 
particle and ^He nucleus previously reported, are shown in table together with 
the estimate of the energy of symmetric nuclear matter (SNM) at the empirical 
equilibrium density po =0.16/m~^, calculated with the FP-AFDMC by simulating 
28 nucleons in a periodic box[48j (see the following sections for more details). 



nucleus 


Eafdmc 


Eafdmc/ A 




4He 


-27.13(10) 


-6.78 


-7.07 


8He 


-23.6(5) 


-2.95 


-3.93 


160 


-100.7(4) 


-6.29 


-7.98 


40 Ca 


-272(2) 


-6.80 


-8.55 


oo 




-12.8(1) 


-16 



Table 4.6: The FP-AFDMC ground-state energy in MeV of ^He, ^He, ^^O and ^OCa 
for the Argonne AV6' interaction [47j. Experimental energies are also reported|130j. 
The estimate of the energy of SNM at equilibrium density, simulated with 28 nucle- 
ons in a periodic box[lS] is also included. 



The AV6' interaction is not sufficient to build the whole experimental binding 
energy[130j of the nuclei considered. This NN interaction gives about 96% of total 
binding energy for alpha particle, 75% for ^He, 79% for ^^O and 79% for ^^Ca. By 
neglecting the missing NN terms, the lack of the TNI seems to reach a saturation 
around the value of ~20%. However the ^^O is unstable to break up into 4 alpha 
particles, and the ^OCa has the same energy of 10 alpha particles. This behavior is 
consistent with the simple pair counting argument of Ref. [131j . 

The estimate of SNM energy as well as the ground-state energy of ^OQa can be 
used to fit the surface energy coefficient in the Weizsacker formula [132j that is given 
as a function of the mass number A and atomic number Z: 

BE{A,Z) _ _ J_ _ Z{Z-1) _ {A - IZf 6_ 
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where Uy gives the volume energy, is the coefficient of the surface energy, ac is the 
term of the Coulomb energy, agym gives the asymmetry energy and 5 describes the 
pairing energy. The AFDMC calculations do not contain the Coulomb term. From 
the comparison of the binding energies per nucleon of symmetric nuclear matter and 
^°Ca, the fitted surface coefficient is 20.5 MeV, not too far from the experimental 
value of 17.23 MeV[l33]. 

4.2 Neutron drops 

A neutron drop is a simple system, useful to benchmark theoretical predictions, and 
significant both for nuclear and astronuclear physics. Several works on neutron drops 
were presented to extract some information for mean field calculations |134l I135| of 
neutron matter. 

In order to test the ability to reproduce the ground-state of neutron drop the 
CP- AFDMC was applied to calculate the energy of a drop as well as the drop 
in both the J=3/2 and J =1/2 states [T35] . and results were compared to the GFMC 
energy. The results indicated that the ground-state of drop was in agreement with 
GFMC within 2%. Despite the agreement on the total energies, the yielded spin-orbit 
splitting (SOS) derived as the difference between the two states with different 
J, was too small compared to the GFMC estimate. This fact was interpreted as a 
problem in dealing with the spin-orbit contribution of the AFDMC algorithm, as 
just suggested by earlier neutron matter calculations (that however will be discussed 
in the following sections). 

In addition, as just discussed in chapter [2] the NN interaction is not sufficient 
to well describe properties of light nuclei and different forms of TNI were proposed 
to build a nonrelativistic Hamiltonian that well reproduces experimental results [63] 
such ground-state energy, density profile, square mean radius and others. The Ur- 
bana IX TNI (UIX) was implemented to well describe properties of nuclei with 
A<4[23j, but it was not satisfactory to reproduce heavier nuclei. An alternate form 
is the TNI Illinois-1 to 5 (ILl to IL5), developed to well describe nuclei up to 
A<8[63], and later employed to compute the ground state energy of nuclei with 
A<10[81J and subsequently up to A=12[35j. However it was recently shown that 
these accurate TNI interactions have some problems in reproducing the properties 
of neutron-rich nuclei like ^'^He. 

CP- AFDMC calculation on neutron matter revealed that Illinois-like (ILx) TNI 
give very unexpected results. In particular, ILx gives a large attractive contribution 
to the energy when the density of neutron matter increases [137 ]. With the same NN 
interaction, at equilibrium density po=0.16 fm~^ all the ILx TNI interactions give 
similar results (within 10%), but when density increases their contribution vary in 
a sensible way. 

Both points are hence discussed. 1) It will be shown that FP-AFDMC better 
reproduces the SOS, which is now comparable to that given by GFMC. 2) In order 
to better clarify the role of the nuclear TNI in pure neutron systems, I will discuss 
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some ground-state energy of neutron drops made of 8 and 20 neutrons at several 
densities. 

The ground-state energy of a neutron drop can be calculated by starting from a 
nonrelativistic Hamiltonian of the form of Eq. 12.11 with the addition of an external 
field: 

^ = - E 1^^' + E + E % + E ^^1^ ' (4-2) 

i i i<j i<j<k 

where the Vij NN interaction is the Argonne AV8' (see section [2T]) . the TNI consid- 
ered were the Urbana IX and Illinois ILl-4 forms (see section [212]) . and the external 
Wood-Saxon well is needed to have a bound state for pure neutron system, otherwise 
unbound. The external field has the following form: 

^-*W = -T^Jb)7^' (4.3) 

where the parameter a =0.65 fm is fixed, while the Vq and R were varied to modify 
the density of the drop. 

The trial wave function has the usual form with the omission of the isospin 
degrees of freedom, and radial orbitals were obtained by solving the Hartree-Fock 
problem with the Skyrme SKM force |138] . The wave functions employed for drops 
with 8 and 20 neutrons are built as for closed shells. For 8 neutrons we only need 
lS'i/2, I-P3/2 and l-Pi/2 orbitals, and with the bigger droplet with 20 neutrons we 
also add ID5/2, 251/2 1^3/2 orbitals. In the case of 7 neutrons, the two states 
with J = 3/2 and J = 1/2 are obtained by dropping the orbitals (/>i,3/2,3/2 and 
01,1/2,1/2 respectively. 

The results of ''n drop with the corresponding SOS contributions for various 
Hamiltonian are reported in table 14.71 The fixed-phase approximation noticeably 
improves the AFDMC results for SOS. In fact, the energy of ^n(l/2"'") is lower 
respect that obtained with the constrained-path [136] and in a better agreement 
with GFMC[63]- The two-body interaction employed is not the same, in fact we use 
the Argonne AV8' while in GFMC the more accurate AV18 is employed. However, 
the SOS given by GFMC with AV8' and UIX is 1.6 MeV, and the overall agreement 
is good. Small differences are present with the IL3 where the SOS is a bit different 
but however well within error bars. 

The first main conclusion is that AFDMC with the fixed-phase can correctly 
build the spin-orbit correlations that seemed to be not satisfactory with the constrained- 
path approximation. 

Let then be examined the ground-state energy as a function of the external well 
parameters of and '^^n drops respectively reported in table 14.81 and 14. 9[ The a 
parameter is fixed to the value of a=0.65 fm, while the Vq and R are varied to change 
the density of the drop. Reported results were obtained using the Argonne AV8' NN 
interaction, and with the addition of various TNI considered as indicated. The TNI 
contribution, obtained by subtracting for each drop the corresponding AV8' result, 
are reported in table 14.101 and I4.11[ 
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TNI 


■''n(l/2+) 


■^n(3/2+) 


AFDMC-SOS 


GFMC-SOS 


UIX 


-33.06(3) 


-31.51(2) 


1.55(5) 


1.5(1) 


ILl 


-35.28(3) 


-32.58(2) 


2.70(5) 


2.8(3) 


IL2 


-35.36(4) 


-32.43(1) 


2.93(5) 


2.8(3) 


IL3 


-36.06(4) 


-32.67(2) 


3.39(6) 


3.6(4) 


IL4 


-35.00(3) 


-32.53(3) 


2.47(6) 


2.4(4) 



Table 4.7: The FP-AFDMC ground-state energy of 7 neutrons in a Wood-Saxon well 
with Vo=-20 MeV, R=3.0 fm and a=0.65 fm. Calculations were performed both with 
J=3/2 and J=l/2 with Argonne AV8' two nucleon interaction and with different 
TNI forces. In the last column the GFMC SOS from Ref. [63] were also reported. 
With the UIX the CP-AFDMC gives a SOS=0.7(2) MeVp6]. All the energies are 
express in MeV. 



Vo 


R 


AV8' 


AV8'/UIX 


AV8'/IL1 


AV87IL2 


AV87IL3 


AV8VIL4 


20 


2.5 


-18.66(4) 


-17.98(4) 


-19.69(4) 


-19.72(4) 


-20.02(4) 


-19.61(4) 


20 


3.0 


-38.51(6) 


-37.55(2) 


-39.76(6) 


-39.72(3) 


-40.16(7) 


-39.59(4) 


20 


3.5 


-57.18(5) 


-56.30(6) 


-58.35(5) 


-58.25(5) 


-58.66(5) 


-58.19(4) 


30 


2.5 


-56.94(5) 


-54.82(5) 


-60.08(4) 


-59.8(1) 


-61.2(2) 


-59.60(6) 


30 


3.0 


-87.96(1) 


-86.11(6) 


-91.2(2) 


-90.98(5) 


-92.15(9) 


-90.69(4) 


30 


3.5 


-115.38(5) 


-113.80(5) 


-117.69(4) 


-117.5(1) 


-118.24(8) 


-117.51(4) 


40 


2.5 


-101.01(2) 


-97.92(3) 


-107.01(8) 


-106.58(9) 


-108.82(9) 


-106.6(1) 


40 


3.0 


-142.81(3) 


-140.11(3) 


-148.07(8) 


-147.7(1) 


-149.5(1) 


-147.52(8) 


40 


3.5 


-177.82(5) 


-175.65(4) 


-181.30(5) 


-181.04(4) 


-182.32(7) 


-180.84(6) 


50 


2.5 


-149.12(7) 


-144.89(7) 


-158.0(3) 


-157.4(1) 


-160.8(2) 


-157.05(9) 


50 


3.0 


-200.46(5) 


-196.77(4) 


-206.95(6) 


-206.81(4) 


-209.3(1) 


-206.42(6) 


50 


3.5 


-242.05(4) 


-239.36(5) 


-246.62(7) 


-246.4(1) 


-248.17(8) 


-246.26(7) 


60 


2.5 


-199.66(7) 


-194.28(7) 


-211.1(2) 


-210.5(2) 


-214.9(3) 


-210.2(2) 


60 


3.0 


-260.14(4) 


-255.59(3) 


-268.9(1) 


-268.5(1) 


-271.7(2) 


-267.92(8) 


60 


3.5 


-308.01(6) 


-304.84(4) 


-313.8(1) 


-313.5(1) 


-315.7(1) 


-313.25(6) 



Table 4.8: AFDMC ground-state energy of 8 neutrons in a Wood-Saxon well as a 
function of parameters Vq (in MeV) and R (in fm) of Vext^ as reported in the first 
two columns. Only the a=0.65 fm parameter was fixed. In the third column AV8' is 
the result of energy with only Argonne AV8' NN interaction, and from the fourth 
to eighth are those of NN plus various TNI forces. All the energies are express in 
MeV. 



The energy of the drop strongly depends on the parameters of external well. The 
first observation is that all Illinois forces give always an attractive contribution to 
the drop while the UIX term is repulsive. This is due to the three pions exchange 
term that only appears in the Illinois forces; this term should be always attractive 
in neutron drops and gives a large absolute contribution to the TNI with respect of 
two pions terms as reported by Pieper et al. in [63j for and neutron drops. 
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R 


AV8' 


AV87UIX 


AV87IL1 
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-202.41(5) 


-208.0(3) 


-202.5(4) 


30 


4.0 


-250.99(5) 


-241.48(7) 


-266.0(2) 


-264.7(1) 


-269.0(2) 


-265.0(2) 


40 
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-221.00(8) 


-203.81(6) 


-258.3(3) 


-255.9(2) 


-265.9(3) 


-256.8(1) 


40 


3.5 


-319.77(7) 


-303.65(9) 


-351.5(3) 


-348.7(1) 


-357.7(1) 


-349.9(3) 


40 


4.0 


-404.51(5) 


-390.96(9) 


-427.7(3) 


-426.0(1) 


-431.9(2) 


-426.0(2) 


50 


3.0 


-341.73(8) 


-317.45(6) 


-398.6(4) 


-396.9(1) 


-408.5(4) 


-397.8(3) 


50 


3.5 


-464.0(1) 


-442.31(3) 


-506.9(3) 


-505.4(4) 


-516.5(3) 


-505.7(3) 


50 


4.0 


-564.93(8) 


-547.7(1) 


-595.4(3) 


-594.1(2) 


-601.9(1) 


-594.2(1) 



Table 4.9: AFDMC ground-state energy of 20 neutrons as explained in table l4^ All 
the energies are express in MeV. 



Vo 


R 


UIX 


ILl 


IL2 


IL3 


IL4 


20 


2.5 


0.68 


-1.03 


-1.06 


-1.36 


-0.95 


20 


3.0 


0.96 


-1.25 


-1.21 


-1.65 


-1.08 


20 


3.5 


0.88 


-1.17 


-1.07 


-1.48 


-1.01 


30 


2.5 


2.12 


-3.14 


-2.86 


-4.26 


-2.66 


30 


3.0 


1.85 


-3.24 


-3.02 


-4.19 


-2.73 


30 


3.5 


1.58 


-2.31 


-2.12 


-2.86 


-2.13 


40 


2.5 


3.09 


-6.00 


-5.57 


-7.81 


-5.59 


40 


3.0 


2.70 


-5.26 


-4.89 


-6.69 


-4.71 


40 


3.5 


2.17 


-3.48 


-3.22 


-4.50 


-3.02 


50 


2.5 


4.23 


-8.88 


-8.28 


-11.68 


-7.93 


50 


3.0 


3.69 


-6.49 


-6.35 


-8.84 


-5.96 


50 


3.5 


2.69 


-4.57 


-4.35 


-6.12 


-4.21 


60 


2.5 


5.38 


-11.44 


-10.84 


-15.24 


-10.54 


60 


3.0 


4.55 


-8.76 


-8.36 


-11.56 


-7.78 


60 


3.5 


3.17 


-5.79 


-5.49 


-7.69 


-5.24 



Table 4.10: The TNI contributions to the drop as a function of parameters Vq (in 
MeV) and R (in fm) of Vext-, as reported in the first two columns. Only the a=0.65 
fm parameter was fixed. The TNI contributes to the total energy were obtained 
by subtracting from the full Hamiltonian with each TNI the AV8' result. All the 
energies are express in MeV. 



Hence it is clear that for pure neutron systems the main physics of TNI comes from 
37r exchange and not from 27r exchange terms. 

A global trend of the Illinois potentials is observed by changing parameters of 
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R 


UIX 


ILl 


IL2 


IL3 


IL4 


20 


3.0 


4.2 


-3.8 


-3.1 


-4.9 


-3.4 


20 


3.5 


4.8 


-7.9 


-6.7 


-9.33 


-7.0 


20 


4.0 


4.9 


-8.0 


-6.91 


-9.1 


-7.0 


30 


3.0 


9.5 


-19.4 


-18.0 


-24.6 


-18.8 


30 


3.5 


10.82 


-19.08 


-17.69 


-23.28 


-17.78 


30 


4.0 


9.51 


-15.01 


-13.71 


-18.01 


-14.01 


40 


3.0 


17.19 


-37.3 


-34.9 


-44.9 


-35.8 


40 


3.5 


16.12 


-31.73 


-28.93 


-37.93 


-30.13 


40 


4.0 


13.55 


-23.19 


-21.49 


-27.39 


-21.49 


50 


3.0 


24.28 


-56.87 


-55.17 


-66.77 


-56.07 


50 


3.5 


21.69 


-42.9 


-41.4 


-52.5 


-41.7 


50 


4.0 


17.23 


-30.47 


-29.17 


-36.97 


-29.27 



Table 4.11: The TNI contributions to the ^^n drop as explained in table HTOl All 
the energies are express in MeV. 



the external well (therefore varying the density of the system). For drop with 
Vo=20 MeV and 30 MeV the larger contribution of TNI is given with R=3.0 fm, 
and in the case of ^0=30 MeV the TNI with R=2.5 fm is similar to that of R=3.0 
fm. When Vq increases to 40 MeV, 50 MeV and 60 MeV the contribution of TNI 
drastically changes and the more bound is given if R is smaller (then when the 
density increase). For "^^n drop with Vo=20 MeV the larger TNI contribution is for 
the higher value of R instead of lower for drop. However when Vq increases to 
30 MeV, 40 MeV and 50 MeV the larger attractive contribution from TNI is given 
when R is smaller. In the case of UIX a similar behavior is observed. In this case 
the TNI contribution is always repulsive, but the absolute value of TNI essentially 
has the same trend of Illinois forces. However, it has to be pointed out that the UIX 
contribution seems to change more slowly than that of Illinois when the density 
of the system increases. These observations suggest that probably the attractive 
contribution given by 37r exchange dependence by density is stronger compared to 
the 27r and central terms. Hence, for smaller densities the repulsion grows as the 
attraction, but when density increases the negative 37r exchange term becomes more 
important than repulsion. This observation is in agreement of observations in pure 
neutron matter calculations [137j . 

The IL2 and IL4 seem to give substantially the same contribution in all the 
systems we simulated, but a small different trend is observed by using a different 
number of neutrons. In fact in the drop the IL2 is as always slightly more at- 
tractive than IL4 but in the case of "^^n the IL4 gives more binding, although the 
absolute values are very close. The IL2 and IL4 have a similar coefficient for 37r 
exchange term, 0.0026 and 0.0021 respectively, and the parameter of the term 
are -0.037 and -0.028. This suggest that probably the contribution of O^'^ is always 
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very small as reported in ref. [63]. 



The ILl is always more bound than IL2, and the difference between them never 
exceeds the 20% of total TNI. The relative difference between ILl and IL2 with 
respect to the total TNI seems to be higher when the Vq of external well is smaller. 
The strength of O^^^ is a bit different but not important as suggested in observing 
the behavior of IL2 and IL4. However IL2 differs from ILl because it has the term 
of that is zero in ILl and this suggest that this term gives always repulsion. 

The IL3 always gives more binding with respect to ILl, IL2 and IL4 and this 
probably is a consequence of the different strength of the ^3^^ parameter that for 
IL3 is 0.0065 respect to 0.0026 of ILl and IL2 and to 0.0021 of IL4. 

The radial density of a given drop can be calculated as a mixed operator |93|: 



where (pQ is the walker distribution in the limit of infinity imaginary time that 
represent the ground-state of the system. Because the trial wave function ■0t does 
not contain any spin correlations, this mixed operator should be not well accurate. 
The mixed-densities for each drop was calculated anyway to show the effect of using 
different Hamiltonians and for different parameters of external well. 

In Fig l4.1l we plot the radial densities of the '^^n drop in an external well with 
Vo=20 MeV, R=3.0 fm and a=0.65 fm for different Hamiltonians. It is clear that 
UIX and IL2 give completely different contributions to the density profile compared 
to the pure two-body interaction. In particular the UIX lowers the density at the 
center of the drop and the IL2 increases it. This refiects the fact that UIX gives 
a repulsive contribution to the drop energy, while IL2 is overall attractive. All the 
Illinois TNI have this peculiarity as seen in Fig l4.2l showing that all the densities 
computed with different Illinois forces are similar. 

The effect of using different TNI with the two-body interaction becomes bigger 
also when the density increases. The Vq parameter of the external well was changed 
to 50 MeV to increase the density at the center of the drop. As it can be seen 
in Fig |4.3l the presence of TNI drastically increases the center density value. The 
density computed with UIX is a bit lower with respect the pure AV8' NN interaction. 
However, the difference with the IL2 case is evident. 

In Fig l4.4l and 14.51 were reported the densities of AV8' plus IL2 and UIX and by 
varying the external weU depth Vo=20 MeV, 30 MeV, 40 MeV and 50 MeV. 

Finally we studied the density varying the R parameter of external well by 
keeping fixed Vo=20 MeV using the AV8'-|-IL2 interaction. As it can be seen in 
Fig l4. 61 there are some non significant changes to the surface density, as expected. 

In order to extract some information about the energy dependence on density 
we define an average density of the drop p as: 



mix 



{iprlfpo) {'4't\iPt) 



(4.4) 



N 



(4.5) 
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0.006 



AV8' 

■ ■ ■ AV8' + UIX 
-- AV8'+IL2 




Figure 4.1: Radial densities of "^^n drop interacting with the AV8' and AV8'+UIX 
and IL2 three-body force in an external well with Vo=20 MeV, R=3.0 fm and a=0.65 
fm. 
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Figure 4.2: The same of Fig |4.1l with AV8' and different Illinois TNI forces. 
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r[fm] 



Figure 4.3: The same of FigOwith Vo=50 MeV. 




Figure 4.4: Radial densities of neutron drop with 20 neutrons interacting with 
AV8'+IL2 for different values of the parameter Vq of the external well with R=3.0 
fm and a=0.65 fm. 
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Figure 4.6: Radial densities of neutron drop with 20 neutrons interacting with 
AV8'+IL2 interaction for different values of the parameter R of the external well 
with Vb=20 MeV and a=0.65 fm. 
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where N is the number of neutrons in the drop, and V=|7ri2Q is the volume of the 
drop, with Rq calculated as 

<rI>= j r^p{r)dr . (4.6) 
The energy as a function of the p for the '^^n is reported in Fig J4.7l 




P [fm '] 



Figure 4.7: Energy as a function of average density p for Hamiltonians with AV8' 
plus IL2 and UIX. Each curve correspond to different values of R parameter of the 
external well. 

A strong dependence of the energy on the parameter R appears. The effect of 
R is related only to the surface of the drop and not to the center, as noted above. 
The UIX and IL2 give completely different contributions for the same R, indicating 
that the effect of different TNI is probably not only due to the different density in 
the center of the drop, but also to the surface effects. 



Chapter 5 



Results: nuclear and neutron matter 



The main challenging and fundamental problem in nuclear physics and nuclear as- 
trophysics is the determination of the equation of state (EOS) of nuclear matter. 
The stability of neutron-rich nuclei jl39] . the dynamics of heavy- ion collisions [6| 1 140] . 
the structure of neutron stars [2], and the simulation of core-collapse supernova |141] . 
all depend on the nuclear matter EOS jl42 ]. 

The calculation of the EOS of nuclear matter models with realistic NN interac- 
tions is considered of primarily importance in astronuclear physics. Starting from 
the famous "Bethe Homework problem" invented by Hans Bethe to test the exist- 
ing many-body theories in dealing with dense and cold hadronic matter [l43t I144j . 
considerable advances have been made p^Sj I146j . However, yet, it is not possible to 
firmly ascertain the degree of accuracy of the approximations one has to introduce 
in any of such theories, and substantial discrepancies still exist among the different 
theoretical estimates of the EOS, the response functions and Green's Functions of 
nuclear matter. 

At present the theoretical uncertainties on the equation of state, coming from 
both the approximations one has to introduce in the many-body methods, and the 
lack of knowledge of the nuclear interaction in the high density regime, do not 
allow for definite conclusions in the comparison with the data from astronomical 
observations |147j . 

In this chapter it will be shown the application of the AFDMC algorithm, within 
the fixed-phase approximation, to calculate the EOS of symmetric nuclear matter 
(SNM) using a simplified NN interaction which still contains full tensor correlations. 
This first work reveals and focuses the problematic differences in principal many- 
body techniques used at present for the EOS calculation and also employed to fit 
the TNI. 

The determination of the pure neutron matter (PNM) EOS obtained with a 
modern and realistic Hamiltonian will be presented, and some conclusion on the 
structure of neutron stars will be given and discussed. 

Finally it will be discussed the application of AFDMC to study the low-density 
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regime of PNM where a superfluid phase is present. 

5.1 Nuclear matter 

As shown in Sec. 14.11 the AFDMC combined with the fixed-phase approximation 
works very weU with nuclear Hamiltonian that contains a tensor force, and permits 
to calculate the ground-state energy of nuclei in good agreement with other accurate 
techniques. 

In order to better understand the status of modern many-body theories and 
their accuracy |148j the AFDMC algorithm can be applied to perform a simulation 
of homogeneous nuclear matter [l8], simply by modifying the trial wave function 
used for the simulation. 

5.1.1 Symmetric nuclear matter, comparison with other techniques 

In order to simulate the nucleonic matter in the normal state, a trial wave function 
using the Fermi-gas antisymmetric part as described in Sec. 13. 71 is employed. As just 
explained, the use of plane waves as single-particle orbitals constrains the number 
of nucleons that can be used in the calculation. The simulation is carried out with 
nucleons in a periodic box the volume of which depends on the density and nucleon's 
number. 

Finite-size effects are then present for two reasons: the first is the dependence 
of the kinetic energy of Fermi-gas to the number of particles, and the second is 
due to the interaction. A Monte Carlo simulation of an infinite system is performed 
in a finite periodic box with size L, and all inter-particle distances are truncated 
to L/2. Usually, tail corrections due to this truncation are estimated with an in- 
tegration of the interaction from L/2 to infinity. However, this is possible only for 
spin-independent terms. In order to include some tail correction to the potential it 
is possible to include the contributions given by neighboring cells to the simulation 
box |137j . Then, each two-body contribution to the potential is given by 

Vp{x,y, 2) = ^ fp (|(x + mL)x + {y + nL)y + {z + oL)z\) (5.1) 

mno 

where L is the box size, and m, n and o are ±1, ±2, ±3... depending on the number 
of shells of the boxes considered. For all the simulations performed, the inclusion 
of 26 additional neighbor (corresponding to m, n and o to be -1, and 1) have 
convergent results. In this way the finite-size corrections due to the interaction should 
be correctly included, but the presence of finite size effects in the kinetic energy 
should be verified. 

The SNM-EOS is calculated by using 28 nucleons (7 for each spin-isospin state) 
interacting with the v'q NN interaction for several densities. The fg interaction is a 
truncation of the Argonne AV8' to the vq level. This interaction was chosen because 
in a recent work different many-body theories using this interaction predicted very 
different results to the SNM-EOS [1^]. 
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Before presenting the AFDMC EOS, let us examine the efficacy of the procedure 
to correct finite-size effects. In order to assess the magnitude of finite-size effects 
different simulations were performed using different numbers of nucleons for the 
highest and lowest densities considered. Results are reported in table 15.11 As it 
can be seen for both p =0.08 fm~^ and p =0.48 fm~^ the results given by 76 and 
108 coincide with that given by 28 nucleons within 3%. Then the effect due to the 
finite-size of the system is relatively small. 





E/A(28) 


E/A(76) 


E/A(108) 


0.5 


-7.64(3) 


-7.7(1) 


-7.45(2) 


3.0 


-10.6(1) 


-10.7(6) 


-10.8(1) 



Table 5.1: AFDMC energies per particle in MeV of 28, 76 and 108 nucleons in a 
periodic box at various densities. 

Then the simulation of 28 nucleons interacting with the fg interaction was re- 
peated for different densities in the range with 0.5<(/9//9o)<3, being po =0.16 fm~^ 
the empirical equilibrium density. The energy per nucleon for considered densities 
obtained with the AFDMC in the fixed-phase approximation are reported in table 
15. 2[ The first observation regard the empirical equilibrium density that is far to 
be reproduced. However this was expected because a pure NN interaction is not 
sufficient and a TNI is necessary to have the correct saturation density in SNM[M]. 



p [fm '■^] 


p/pQ 


E/A 


0.08 


0.5 


-7.64(3) 


0.12 


0.75 


-9.81(4) 


0.16 


1.0 


-11.5(1) 


0.20 


1.25 


-13.0(1) 


0.24 


1.5 


-13.73(7) 


0.28 


1.75 


-14.1(2) 


0.32 


2.0 


-14.0(3) 


0.36 


2.25 


-13.5(3) 


0.40 


2.5 


-12.7(2) 


0.44 


2.75 


-11.7(2) 


0.48 


3.0 


-10.6(1) 



Table 5.2: AFDMC energies per particle in MeV of 28 nucleons in a periodic box at 
various densities. 

AFDMC results are then compared with two of mainly used many-body theories |48j. 
namely the Brueckner-Hartree-Fock (BHF) and the Fermi Hyper Netted Chain in 
the Single Operator Chain (FHNC/SOC). 
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The BHF method essentially recasts the Hamiltonian to a single-particle operator 
and the residual part is evaluated by solving the Brueckner-Bethe-Goldstone (BEG) 
equation. The accuracy of the results obtained with BBC depends on the expansion 
level adopted; the perturbative expansion of BBC is not convergent, but the cluster 
diagrams can be grouped with the number of independent hole lines |1U8] . Several 
attempt to study the convergence of hole- line expansion were recently explored |1 49 1 

USD]. 

The Correlated Basis Functions (CBF) theory is particularly well suited to deal 
with strongly interacting systems, that are described with a many-body wave func- 
tion written with some correlation operators acting on the non-interacting state. 
Unfortunately the energy evaluation of the CBF state is not easy to do. The Vari- 
ational Monte Carlo (VMC) could be used to numerically solve the integration of 
CBF state, but the complexity of wave function make impractical the applicability 
of VMC when the number of nucleons increases. Another possibility to calculate 
the CBF energy is given by the FHNC/SOC theorv [T08] . Unfortunately the solu- 
tion of FHNC/SOC equations provides the summation of an infinite class of cluster 
diagrams like that called " elementary diagrams" . Hence the CBF is a variational the- 
ory that can be solved with the FHNC/SOC that however is not exact because the 
summation is not complete and some elementary diagrams are not fully considered. 

The AFDMC EOS is shown in Fig. 15.11 where different EOS obtained with 
FHNC/SOC and BHF are also reported[ll8j. 

By assuming that FP- AFDMC for SNM calculation is as accurate as shown for 
nuclei, and that finite size-effects are small as suggested by table 15.11 the compar- 
ison of the various EOS suggests the following comments: FHNC/SOC leads to an 
overbinding at high density. A similar indication was found by Moroni et al. [151j 
after a DMC calculation of the EOS of normal liquid ^He at zero temperature, 
with a guiding function including triplet and backflow correlations. The comparison 
with the equivalent FHNC/SOC calculations of Refs. [152\ 1153] have shown similar 
discrepancies. 

On the other hand there are two main intrinsic approximations in variational 
FHNC/SOC calculations, which violate the variational principle. The first one con- 
sists in neglecting a whole class of cluster diagrams, the elementary diagrams, which 
are not summed up in the FHNC integral equations. The calculation of the lowest 
order diagram of this class shows a substantial effect from this diagram represented 
by the blue arrows in Fig. 15. H and bring the FHNC/SOC estimates very close to 
AFDMC results. An accurate study of corrections given by higher diagrams was 
performed by Wiringa |154j showing that their effects are not negligible. The second 
approximation is related to the non-commutativity of the correlation operators en- 
tering the variational wave function. The only class of cluster diagrams contributing 
to such non-commuting terms, which can be realistically calculated, is that charac- 
terized by single operator chains (SOC approximation). It is believed that such an 
approximation is reliable in nuclear matter, but there is no clear proof of this. 

BHF calculations of Ref. [148j includes corrections up to the two hole-line expan- 
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Figure 5.1: Equation of state of symmetric nuclear matter calculated with differ- 
ent methods. Red circles represent AFDMC results with statistical error bars and 
the green line is the fitted functional. Dashed lines correspond to calculations per- 
formed with other methods [148] (blue line with squares: FHNC/SOC; magenta with 
diamonds: BHF). Blue triangles represent the FHNC/SOC energies corrected by in- 
cluding the low order of elementary diagrams. Blue arrows show the corresponding 
energy shift, which increases at higher densities. Figure taken from Ref. |48j . 



sion and predict an EOS with a shallower binding than the AFDMC one. It has been 
shown for symmetric nuclear matter, using the Argonne AVIS and AV14 potentials, 
that contributions from three hole-line diagrams add a repulsive contribution up to 
~ 3MeV at densities below /?n |149j . and decrease the energy at high densities [155]. 
Such corrections, if computed with the v'q based on AV8' potential, would probably 
preserve the same general behavior, and bring the BHF EOS closer to the AFDMC 
one. Therefore, the AFDMC calculations show that the two hole-line approximation 
used in Ref. |148j is too poor, particularly at high density. 

Although the EOS calculated with the AFDMC started from a simplified Hamil- 
tonian, it is particularly important because it emphasizes some limitations of the 
FHNC/SOC calculation that was used to fit the TNI on nuclear matter. At the 
present the TNI forces are fitted using the GFMC to simulate light-nuclei, in the 
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addition to the FHNC/SOC calculation to reproduce the equilibrium density of 
nuclear matter; differences shown between AFDMC and FHNC/SOC open some 
questions about the physical meaning of such potentials in the bulk matter calcula- 
tions, and stress the fact that more work should be addressed to study the role of 
TNI, especially at higher densities. 

The computed EOS with AFDMC is fitted with a functional form as following: 

^ = ^ + a{~p-Pof + (i{~p-Pof, (5.2) 

where p is the nucleon density respect to the equilibrium density 

P=-, (5.3) 
Po 

and the Eq, a, {3 and pq coefficients are given by the fit: 



-14.04(4)Mel^ 



A 

a = 2,.m{<d)MeV , 

13 = -0.44(8)Mey, 

po = 1.83(1). (5.4) 



The nuclear matter compressibility K is given by 

'diE/A) 



K = 9pl 



dp- 



n2 



(5.5) 



po 



The value obtained with the functional form of Eq. 15.21 is K~190 MeV. The 
fitted EOS obtained with the AFDMC could be used to calibrate other many-body 
theories. 

The FP-AFDMC is an accurate method to calculate the EOS of SNM, but 
it could be easily generalized to microscopically simulate the asymmetric nuclear 
matter particular important to study neutron stars. At the present the asymmetric 
nuclear matter can be studied using mean field calculations [1561 [157]. 
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5.1.2 Equation of state of symmetric nuclear matter with the Ar- 
gonne AV6' interaction 

The AFDMC with the fixed-phase approximation was used to calculate the equation 
of state of nuclear matter using the realistic Argonne AV6' NN interaction. 

The Argonne AV6' is obtained by dropping the spin-orbit operators from the 
AV8'. The spin-orbit does not contribute to the S wave in the NN scattering as well 
as the P wave in the spin-singlet ^Pi channel, but it is important in differentiating 
between ^-Po,i,2 channels. Then the AV6' was obtained to preserve the deuteron 
energy by including the energy difference from the AV8' in the central part in the 
5, T=1,0 channel. 



p [fm 


p/po 




E/A(AV6') 


0.04 


0.25 




-5.01(4) 


0.08 


0.5 


-7.64(3) 


-8.14(4) 


0.12 


0.75 


-9.81(4) 


-10.86(3) 


0.16 


1.0 


-11.5(1) 


-12.77(7) 


0.20 


1.25 


-13.0(1) 


-14.67(4) 


0.24 


1.5 


-13.73(7) 


-16.0(2) 


0.28 


1.75 


-14.1(2) 


-16.4(1) 


0.30 


1.875 




-16.6(1) 


0.32 


2.0 


-14.0(3) 




0.36 


2.25 


-13.5(3) 


-16.9(2) 


0.40 


2.5 


-12.7(2) 


-16.7(1) 


0.44 


2.75 


-11.7(2) 


-15.9(1) 


0.48 


3.0 


-10.6(1) 


-14.3(2) 


0.52 


3.25 




-13.0(2) 



Table 5.3: AFDMC energies per particle in MeV of 28 nucleons in a periodic box at 
various densities. The reported ground state energy were obtained by considering the 
Argonne AV6' interaction and the AV8' without the spin-orbit contribution (called 

The AFDMC was employed to simulate 28 nucleons interacting with the AV6' 
as described in the previous section. The tail corrections to the potential were also 
included. The energy of 28 nucleons in a periodic box are reported in table 15.3] were 
they are also included the energies of Ref. [48] obtained with the v'q interaction. 

It is interesting to note how the reprojection of the AV8' to the AV6' by dropping 
the spin-orbit significantly changes the EOS of the symmetric nuclear matter as it 
can also be seen in Fig. 15. 2i The spin-orbit term of AV8' is very small in the alpha- 
particle, and its contribution is -4.75 MeV with respect to the total potential energy 
of -128.25 MeVpS]. The alpha-particle total energy is -25.14 MeV with the AV8' 
and -26.15 MeV with the AV6'[58j. and the energy with the AV8' without the spin- 
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0.48 



Figure 5.2: The AFDMC equation of state of symmetric nuclear matter. The energies 
were evaluated by simulating 28 nucleons in a periodic box interacting with different 
Hamiltonians. The red points (taken from Ref. [48]) were obtained with the Argonne 
AV8' NN interaction without the spin-orbit contribution (called v'q), while the black 
points were calculated with the Argonne AV6' interaction. The lines correspond to 
the fitted functionals. 



orbit is -21.18 MeV^25j. Then the difference between the v'q and the AV6' is of about 
1.2 MeV per nucleon in the alpha-particle ground state. In the nuclear matter at 
densities p =0.08 fm~^, 0.12 fm~^ and 0.16 fm~^ the difference between the v'q and 
the AV6' energies is respectively of 0.5 MeV, 1.05 MeV and 1.27 MeV, not far from 
that in the alpha-particle. When the density increases the spin-orbit contribution 
becomes larger, and also the saturation density increases from 0.29 fm-^ to 0.34 
fm~^. However this effect could also be due to the different strength of the central 
term in the S, T=10 channel between the two interactions. For densities higher than 
~0.36 fm~^ the behavior of the two EOS seems to be similar with a constant shift. 



The calculated EOS with AV6' is fitted with the same functional form of Eq. 15.21 
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as in the previous section: and in this case the Eq, a, (3 and pQ coefficients are: 

^ = -17.0(l)Mey, 

a = 3.1(l)Mey, 
/3 = -0.22(7)Mey, 

po = 2.09(2). (5.6) 

The nuclear matter compressibihty for the AV6' results to be K~244 MeV with the 
functional form of Eq. | 



5.1.3 Equation of state of isospin-asymmetric nuclear matter with 
the Argonne AV6' interaction 

In order to test the reliability of AFDMC to deal with the asymmetric nuclear matter 
calculation, several simulations with N ^ Z were performed. 

In according to the mass formula of Eq. 14.11 the asymmetry parameter a can be 
defined as 

and, in order to have a closed-shell configuration in the trial wave function, only 
particularly (A^, Z) pairs are admitted. The trial wave function has the same struc- 
ture of that used for symmetric nuclear matter, and both A^ and Z must be numbers 
corresponding to a closed shell of plane waves. 

The finite-size effects for the potential were included as described in the symmet- 
ric nuclear matter section, and for density /3=0.08 fm~^, 0.16 fm~^ and 0.24 fm~^'^ 
the energy dependence to the a parameter were studied. 

The AFDMC energies of asymmetric nuclear matter for the considered densities 
are reported in table 15.41 and plotted in Fig. 15.31 As it can be seen the energy 
dependence to the a parameter becomes stronger when the density increases. In 
particular for a > 0.5 there is a strange trend of the energy, probably due to the 
finite-size effect in the kinetic energy. The potential is corrected to avoid these effects, 
but the kinetic energy could still suffer of the fact that different shells are considered. 
This trend is less evident for small values of a probably because in the region of 
a < 0.5 the energy dependence to the asymmetry parameter is weaker. 

The energy per nucleon as a function of a is plotted if Fig. 15.41 where for each 
density is also displayed the quadratic functional of E(a) fitted on the SNM (a = 0) 
and PNM (a = 1) points. The fit at the equilibrium density p=0.1G fm~^ represent 
the coefficient of the asymmetry energy in the Weizsacker formula of Eq. 14. 1[ 

The energy dependence to the asymmetry parameter in the infinite limit is 



2 



E{A,Z) _ {A - 2Z) 

— flu + O-syra 

and by taking -o^=E/N(a=0)=-12.77 MeV, the fit gives asym=27.33 MeV while 
the experimental value is 23.2875 MeV |133] . 
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N 


P 


a 


E/A(p=0.08fm-^) 


E/A(p=0.16fm-^) 


E/A(p=0.24fm-^) 


14 


14 


0.0 


-8.14(4) 


-12.77(7) 


-16.0(2) 


66 


38 


0.27 


-6.26(8) 


-10.29(3) 


-12.3(1) 


38 


14 


0.46 


-4.35(6) 


-7.07(7) 


-7.8(2) 


66 


14 


0.65 


-0.21(4) 


-0.36(6) 


0.4(4) 


14 


2 


0.75 


0.98(2) 


0.99(6) 


2.25(7) 


38 


2 


0.9 




8.36(2) 




66 


2 


0.94 




12.22(2) 




14 





1.0 


9.43(1) 


14.56(2) 


20.23(3) 



Table 5.4: AFDMC energies per particle in MeV of asymmetric nuclear matter in 
a periodic box at various densities. The number of neutrons and protons (A^ and 
Z) and the corresponding asymmetry parameter a is indicated. The Hamiltonian 
contains the Argonne AV6' NN interaction. 




0.2 0.4 0.6 0.8 1 



a 

Figure 5.3: The AFDMC energy of asymmetric nuclear matter as a function of the 
asymmetry parameter a at different densities as indicated. The considered Hamil- 
tonian is the AV6', and the estimate of the energies were obtained by simulating 
different set of (A^, Z) in a periodic box. 
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Figure 5.4: The AFDMC energy of asymmetric nuclear matter at different densities 
as a function of the asymmetry parameter. 



The effect of the tensor term is essential in the asymmetry parameter calcula- 
tion. In fact a previously AFDMC using a f4-like interaction predicted 0^^^=36.4 
MeV |118] that is sensibly larger with respect to both the experimental fit and to 
the AFDMC calculation with the AV6' interaction. 

However, an accurate calculation of the asymmetric nuclear matter equation of 
state requires a different treatment of the finite-size effects due to the kinetic energy. 
One possible way to make shell effects milder is the use of the twisted averaged 
boundary conditions [llOj . 
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5.2 Neutron matter in the normal phase 

An approximation of matter as it is found in neutron stars is a system made of 
neutrons only, called pure neutron matter (PNM) . The AFDMC algorithm combined 
with the constrained-path approximation was already used by Sarsa et al. to study 
the equation of state of PNM at zero temperature [137J. The Hamiltonian considered 
contains both a realistic AV8' NN and the Urbana IX TNI interactions; at present 
such Hamiltonian is the most commonly used to calculate properties of both SNM 
and PNM. 

The CP- AFDMC proved to give a very satisfactory result for PNM calculations 
with a NN and TNI interaction, but some problem were encountered in the eval- 
uation of the spin-orbit contribution as just discussed in the neutron drop section 
I4.2[ The inclusion of spin-backflow correlations reduced the discrepancies |158J, but 
a detailed study just considering a pure NN interaction emphasized this problem 
of AFDMC in dealing with the spin-orbit correlations [159]. A similar behavior was 
found by comparing the AFDMC |137j with the GFMC evaluation of the energy of 
14 neutrons in a periodic box [36] I37j . By considering the same Hamiltonian with 
the same box truncation used in GFMC calculation, the AFDMC overestimates the 
energy of 14 neutrons with a AV8' interaction. 

The FP- AFDMC overcomes the problem of the spin-orbit correlations as shown 
in table 15. 5i Without tail corrections the CP- AFDMC energy of 14 neutrons at 
p =0.16 fm-3 is 20.32(6) MeV compared to 17.00(27) MeV given by UC-GFMC, 
while the FP-AFDMC energy is 17.586(6) MeV, in better agreement (within 3%) 
with the UC-GFMC. For higher densities reported in the table it should be noted 
that the CP-GFMC significantly differs from UC-GFMC, probably because the 
Fermion sign problem in that case is more severe [36]. and the convergence of the en- 
ergy can be very problematic, because the evolution of the energy can be carried out 
only for a very short imaginary-time. These shortcomings could introduce some spu- 
rious effects limiting the accuracy of GFMC for the neutron matter calculation [117j . 
Anyway this comparison in addition to the SOS calculation in neutron drops demon- 
strates the high accuracy of FP-AFDMC to deal with neutron systems interacting 
with realistic NN Hamiltonians. 

The FP-AFDMC was applied to study the PNM by simulating periodic boxes 
containing a different number of neutrons interacting with the Argonne AV8' NN 
potential, including the finite-size corrections to the potential energy as described 
in the above section of symmetric nuclear matter. All FP-AFDMC results, as well 
as those of Ref. [T37] with CP- AFDMC are reported in table El Where available, 
the CP- AFDMC calculation using the same Hamiltonian and the same corrections 
are included for comparison. 

The first consideration concerns the inclusion of finite-size effects that sensibly 
reduce the gap between the FP-AFDMC and CP-AFDMC with 14 neutrons. In fact 
without this corrections for AV8' at p =0.16 fm~'^ the difference is less than 14%, 
while in this case is less than 3%. For higher densities the discrepancy never exceed 
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p [fm '■^] 


FP-AFDMC 


CP-AFDMC 


CP-GFMC 


UC-GFMC 


0.04 


6.69(2) 




6.43(01) 


6.32(03) 


0.08 


10.050(8) 




10.02(02) 


9.591(06) 


0.16 


17.586(6) 


20.32(6) 


18.54(04) 


17.00(27) 


0.24 


26.650(9) 




30.04(04) 


28.35(50) 



Table 5.5: FP-AFDMC energies per particle of 14 neutrons interacting with the 
AV8' NN interaction in a periodic box at various densities. The CP-AFDMC of Ref. 
[137], the constrained path CP-GFMC and the unconstrained UC-GFMC of Ref. 
|36j are also reported for a comparison. All the energies are expressed in MeV. 



p [fm^l 


FP-AFDMC(14) 


FP-AFDMC(38) 


FP-AFDMC(66) 


CP-AFDMC(14) 


0.12 


11.843(9) 


10.69(2) 


11.99(3) 


12.32(5) 


0.16 


14.566(7) 


13.23(3) 


15.02(7) 


14.98(7) 


0.20 


16.76(1) 


15.13(3) 


17.11(4) 


17.65(7) 


0.24 


20.05(2) 


17.66(3) 


19.87(4) 




0.28 


23.142(8) 


20.76(2) 


23.36(4) 




0.32 


26.038(8) 


23.88(3) 


26.74(4) 


27.3(1) 


0.36 


30.28(2) 


27.35(2) 


30.20(5) 




0.40 


33.344(7) 


30.87(6) 


34.32(4) 


35.3(1) 


0.48 






43.14(5) 




0.56 






53.34(6) 




0.64 






64.90(6) 




0.80 






91.28(3) 





Table 5.6: FP-AFDMC energies per particle of 14, 38 and 66 neutrons interacting 
with the AV8' NN interaction in a periodic box at various densities. The finite- 
size effects due to the NN truncation are included. The CP-AFDMC energy of 14 
neutrons with the same interaction and finite-size corrections [137] are also reported 
for a comparison. All the energies are expressed in MeV. 

the 5% respect to the total energy. 

Some finite-size effects are still present, as it can be seen from the energies ob- 
tained by varying the number of neutrons. The same behavior is consistently followed 
at each density, and typically £'(38) < i?(14) < £'(66), although small deviations are 
present, and in some case £(14) > £(66). This trend directly follows from behavior 
of the kinetic energy of the Fermi gas. As discussed above, finite-size effects due to 
the kinetic energy could however be included using the twisted averaged boundary 
conditions [llOj. It should also be noted that the finite-size contributions due to the 
spin-orbit potential are neglected, and the spin-orbit potential is obtained by a sum- 
mation of pairs inside the simulation box. The effect of this approximation should 
be small, but remains to be studied. 

For PNM the Urbana UIX TNI force reduces to a NN spin-independent interac- 
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p [fm ■^J 


FP-AFDMC(14) 


CP-AFDMC(14) 


JSB-CP-AFDMC(14) 


0.12 


14.12(1) 


14.80(9) 




0.16 


18.76(1) 


19.76(6) 




0.20 


23.57(2) 


25.23(8) 




0.32 


45.51(1) 


48.4(1) 


46.8(1) 


0.40 


65.86(2) 


70.3(2) 




P 


FP-AFDMC(38) 


CP-AFDMC(38) 




0.12 


13.09(1) 


13.96(5) 




0.16 


17.84(2) 


18.67(6) 




0.20 


22.76(4) 


24.7(1) 




0.32 


46.84(3) 


46.8(2) 




0.40 


70.07(5) 


76.3(2) 




P 


FP-AFDMC(66) 


CP-AFDMC(66) 


JSB-CP-AFDMC(66) 


0.12 


14.36(1) 


15.26(5) 




0.16 


19.64(1) 


20.23(9) 




0.20 


24.71(3) 


27.1(1) 




0.32 


49.94(3) 


54.4(6) 


52.9(2) 


0.40 


74.2(1) 


81.4(3) 





Table 5.7: FP-AFDMC energies per particle of 14, 38 and 66 neutrons interacting 
with the AV8'+UIX interaction in a periodic box at various densities compared 
with the available CP-AFDMC ones of Ref. ^137]. The available CP-AFDMC using 
spin-backflow correlations taken from Ref. |158j are also reported. All the energies 
are expressed in MeV. 

tion as explained in Sec. 13.51 Moreover it can be easily included into the propagator. 
The finite-size corrections due to the UIX are not included because they increase the 
computational time and rather cumbersome. The FP-AFDMC results compared to 
available CP-AFDMC of Ref. |137j are reported in table 15. results were obtained 
by simulating different systems with 14, 38 and 66 neutrons. For 14 and 66 neutrons 
with the same Hamiltonian the CP-AFDMC using spin-backflow correlations are 
also reported. The effect of using the flxed-phase approximation rather then the 
constrained-path to the spin-orbit evaluation are once more evident, also comparing 
the FP-AFDMC with the available CP-AFDMC using spin-backflow correlations. 

To better show the differences between the FP-AFDMC and CP-AFDMC some 
results are plotted in Fig. 15.51 

All the FP-AFDMC results for 14, 38, 66 and 114 neutrons interacting with the 
AV8'-|-UIX Hamiltonian are summarized in table 15.81 where it were also reported 
the extrapolated CP-AFDMC energies to the infinite system. 

The obtained results suggest some consideration: finite-size effects are important, 
especially using 14 or 38 neutrons. However, in the range of considered densities the 
energy of 66 neutrons is always very close to that of 114 within few percent, indicat- 
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Figure 5.5: The FP-AFDMC and CP-AFDMC [HZ] equation of state for AV8' and 
AV8'+UIX. See the legend for details. The EOS of CP-AFDMC il37j and FP- 
AFDMC with AV8' were calculated by simulating 14 neutrons in a periodic box, 
while in the case of AV8'+UIX 66 neutrons were used. 





FP-AFDMC 
(14) 


FP-AFDMC 
(38) 


FP-AFDMC 
(66) 


FP-AFDMC 
(114) 


CP-AFDMC 

(oo) 


0.12 


14.12(1) 


13.09(1) 


14.36(1) 


15.13(1) 


15.5 


0.16 


18.76(1) 


17.84(2) 


19.64(1) 


20.73(2) 


20.6 


0.20 


23.57(2) 


22.76(4) 


24.71(3) 


25.75(6) 


27.6 


0.24 


30.51(1) 


29.24(3) 


31.46(3) 


32.85(5) 




0.28 


37.91(1) 


37.46(1) 


40.19(3) 


41.83(4) 




0.32 


45.51(1) 


46.84(3) 


49.94(3) 


51.83(1) 


55.6 


0.36 


56.336(9) 


57.84(2) 


63.5(1) 


63.5(2) 




0.40 


65.86(2) 


70.07(5) 


74.2(1) 


76.74(6) 


83.5 


0.48 






105.9(1) 






0.56 






145.3(2) 






0.64 






193.5(1) 






0.80 






313.7(2) 







Table 5.8: FP-AFDMC energies per particle of 14, 38, 66 and 114 neutrons inter- 
acting with the AV8'-|-UIX interaction in a periodic box at various densities. The 
extrapolated CP-AFDMC energy for oo-neutrons of Ref. [137j are also reported. All 
the energies are expressed in MeV. 
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ing a possible saturation near the infinite limit. This behavior was observed after the 
study of finite-size effects using the periodic box FHNC (PB-FHNC) technique [160], 
and applied to extrapolate to the infinite limit the CP-AFDMC results reported in 
table I5.8[ The PB-FHNC might also be used to correct the truncation of inter- 
particle distances entering in the spin-orbit and TNI potential; anyway the deep 
analysis of Sarsa et al. |137j suggests that the energy of the system in the infinite 
limit is somewhere between the energies of 66 and 114 neutrons. 

In Fig. 15.61 we plot the FP-AFDMC equation of state obtained as the energy of 
66 neutrons, and the variational calculation of Akmal et al. of Ref. [85j, where the 
AV18 NN interaction was used in combination with the Urbana UIX TNI. As it can 
be seen both AV8' and AV18 essentially give an EOS with the same behavior, but 
the addition of the TNI adds some differences, in particular at higher densities. 




0,2 0,4 0,6 0,8 



Figure 5.6: The FP-AFDMC equation of state evaluated by simulating 66 neutrons 
in a periodic box; the AV8' and AV8'+UIX Hamiltonians were considered. The 
EOS were compared with the variational calculations of Ref. [85j using the AV18 
and AV18+UIX Hamiltonians. See the legend for details. 



5.2.1 Neutron star structure 

The PNM results reported in table 15.81 were used to predict some property of a 
neutron star. The structure of a non-rotating neutron star can be fully determined 
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by solving the Tolman-Oppenheimer-VolkofF (TOV) equations [U [21 11]: 

dP _ G[m{r) + ATTr^P{r) /c^] [e(r) + P(r) /c^] 
dr r[r — 2Gm{r) / c^] ' 

= 47r6(r)r2 , (5.10) 

where P{r) and e(r) are the pressure and the energy density of the matter, m{r) 
is the gravitational mass enclosed within a radius r, and G is the Gravitational 
constant. To solve the TOV, the EOS of baryonic matter is needed, whose main 
bulk is given by the neutron matter physics. 

It is commonly accepted that in the lower density regime where p ~ po the 
matter inside a neutron star is dominated by nucleons and leptons, but when den- 
sity increases other hyperonic species may play an important and not negligible 
role in the physics of the neutron star. The inclusion of hyperons in the EOS is 
far to be accessible because the of lack of knowledge about the nucleon-hyperon 
and hyperon-hyperon interactions; only few experimental data were extracted from 
nucleon-hyperon scattering and some potential were proposed [IJJ. 

As a first approximation, the PNM equation of state calculated with 66 neutrons 
with the FP-AFDMC was employed to solve the TOV equations. For this purpose 
the results of table 15.81 were fitted to the functional form 

E{p) = ap" + (3p^ , (5.11) 

where E{p) is the energy per neutron as a function of the density, and a, a, /3 and 
h are parameters given by the fit. 

Using the fitted energy density e(p) defined as 

e{p) = pE{p) (5.12) 

and the pressure P{p) at zero temperature given by 

P(rt=P^^ (5.13) 

it is possible to solve the TOV equation using as the initial condition the pc pa- 
rameter that represents the density at the center of the star. The distance from the 
center at which P = is the radius R of the star. From the value of the radius it is 
possible to calculate the gravitational mass M defined as 

rR 

M{R) = / dr47re(r)r2 . (5.14) 

The gravitational mass in unit of solar masses Mq as a function of pc or R can 
be directly compared to the experimental observations of neutron stars[3]. The TOV 
solutions are plotted in Fig. 15.71 and 15.81 together with the TOV solutions given by 
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Figure 5.7: Neutron star gravitational mass as a function of the radius given by the 
AFDMC and VCS taken from Ref. [85]. The considered Hamihonians are the AV8' 
and AV8'+UIX for AFDMC, and AV18 and AV18+UIX for VCS. 



the variational chain summation (VCS) of Akmal et al. from Ref. [85], where the 
full AV18 as NN interaction instead of the AV8' was considered. 

In the M-R diagram it is possible to understand the dependence of the gravi- 
tational mass by the radius of the star. The predicted maximum mass and radius 
are not too far from observations despite the great simplifications introduced, and 
the absence of important physics coming from the inclusion of hyperons and quark 
matter. The global behaviour is similar for both the AFDMC and VCS, especially 
in the case of the full Hamiltonian; in particular, the AFDMC predicts a maximum 
mass and radius slightly higher with respect to the VCS with the addition of TNI, 
while the maximum mass is lower for pure NN Hamiltonian. 

In the M-pc diagram the dependence to the central density pc of the gravitational 
mass of the star is displayed, pc is a parameter used to solve the TOV equations but 
can be determined by measurements of the maximum mass. Also in this case the 
AFDMC and VCS results are slightly different, especially in the case of pure NN 
interaction. 

The predicted maximum gravitational masses, with the corresponding central 
density and radius obtained by the TOV solution are shown in table 15.91 The results 
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Figure 5.8: Neutron star gravitational mass as a function of the central density 
Pc calculated with the AFDMC and VCS taken from Ref. [85j . The considered 
Hamiltonians are the same of Fig. 15. 7[ 



of Akmal et al. taken from Ref. [85j are also included for sake of comparison. The 
properties computed in this work are not so far with respect to VCS calculations. 
The different methods considered give almost the same gravitational mass if the 
only NN interactions is used, and the difference does not change much after the 
addition of the UIX interaction. 

It has again to be stressed the fact that a better knowledge and investigation 
of the TNI structure as well as an accurate method to solve the Hamiltonian is of 
primary importance to predict neutron star properties. 
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Hamiltonian 


method 




Pc [fm^^] 


R [Km] 


AV8' 


AFDMC 


1.64 


1.89 


7.9 


AV8'+UIX 


AFDMC 


2.45 


0.91 


11.1 


AV18 


vcs 


L68 






AV18+UIX 


vcs 


2.39 






AV18+(5v+UIX* 


vcs 


2.21 


1.14 





Table 5.9: Maximum gravitational masses, in Mq, the corresponding central density 
Pc and the radius. The AFDMC results are obtained from the EOS of 66 neutrons 
using the indicated Hamiltonians, while the VCS results are taken from Ref. |85j . 

5.2.2 Pure neutron matter and Illinois TNI forces 

The simulation of PNM using the Illinois TNI forces [63] was performed. There are 
signs of some problems [137] similar to those discussed in the neutron drop section 
above. The 2>tt exchange term seems to overbind the neutrons when density increases, 
and the phenomenological central term is probably not repulsive enough. 

The Illinois forces were fitted using the GFMC for light nuclei; the IL2 was 
extended up to ^^C and provided binding energies in very good agreement with 
experimental data[35]. Then, in addition to the AV8' interaction, the IL2 was con- 
sidered, and by means of the FP-AFDMC several simulations of 66 neutrons as 
described in the above section were performed by varying the strength of Vstt and 
Vr to see the effect of the behavior of the PNM-EOS. Results are reported in table 
[5lO] and [5Tn 
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Table 5.10: FP-AFDMC energies per particle calculated with 66 neutrons and 
AV8'-|-IL2 interaction in a periodic box at various densities. The strength of the 
Vstt was varied from the original IL2 also reported. All the values are expressed in 
MeV. 

As it can be seen in Fig. 15. 9^ the effect of the 37r exchange term of IL2 is dominant 
and by varying its strength the EOS drastically changes. Only if < 0.01 MeV 
the EOS is comparable to that obtained using the UIX TNI. 

The phenomenological An part also plays a very important role, but the repulsion 
given by this term is not sufficient to contrast the attraction given by the 37r exchange 
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Table 5.11: FP-AFDMC energies per particle calculated with 66 neutrons and 
AV8'+IL2 interaction in a periodic box at various densities. The strength of the 
Vr was varied from the original IL2 that is ^^^=0. 00705 MeV. All the values are 
expressed in MeV. 
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Figure 5.9: AFDMC equation of state calculated with 66 neutrons interacting with 
AV8'+IL2 interaction. Different EOS were obtained by varying the strength of the 
Astt parameter of IL2. The AV8'+UIX is also reported for a comparison. The A^.,^ 
values are expressed in MeV. 
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term. As it can be seen in Fig. lS.lOl the Afi parameter was increased from the original 
IL2 value of 0.00705 MeV up to 0.0105 MeV, but the EOS is completely different, 
and still far to be comparable with the equation of state computed with the UIX 
instead. 
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Figure 5.10: AFDMC equation of state calculated with 66 neutrons interacting with 
AV8'+IL2 interaction. Different EOS were obtained by varying the strength of the 
Aji parameter of IL2. The AV8'+UIX is also reported for a comparison. The A^t^ 
values are expressed in MeV. 



5.3 Superfluidity in low-density neutron matter 

It has been well established that neutron matter is in a superfluid phase in the 
low-density regime in the region of the crust and in the inner part of a neutron star. 

In this low-density region neutron superfluidity is expected to be mainly due 
to the attractive ^Sq channel, as suggested by the trend of NN phase shifts in 
each scattering channel. The occurrence of the well known virtual neutron-neutron 
channel determines the possibility of a formation of a pairing condensate at low 
density p.61j . The physics of ^So neutron superfluid is important to be understood 
in order to understand several phenomena that occur in the inner crust of a neutron 
star which may in turn be related to the formation of glitches [1621 [163]. 

In recent years several works were performed to estimate the BCS gap for modern 
NN interaction. Results essentially confirm these expectations ^164* il65j . Several 
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groups are in a close agreement on the gap values, converging to an expectation of 
about 3 MeV at a Fermi momentum close to kp ~0.8 fm~^ [1661 11671 1168j . All these 
calculations consider a bare NN interaction as the pairing force. It has been pointed 
out that the screening by the medium of the interaction could strongly reduce the 
pairing strength in this channel |169[ \1&2\ \nO\ I171| . The many-body calculation of 
the pairing gap exactly including all the correlations given by the NN interactions 
is of primary relevance, and the AFDMC can be used to accurately perform a 
calculation of the gap[50|. 

The AFDMC can be used to describe a superfluid system using a BCS trial wave 
function as described in Sec. 13.71 The BCS coefficients appearing in the Pfaffian 
of paired orbital wave functions can be taken from the BCS amplitudes from a 
correlated basis functions (CBF) calculation [50 j. For even N number, all neutrons 
are paired in the Pfaffian, whereas for odd N values the configuration of the unpaired 
neutron providing the best energy must be found. 

The pairing gap energy for odd N values is defined as: 

A(iV) = E{N) - ^[E{N - 1) + ^(iV + 1)] , (5.15) 

and several calculations using different values of N must be performed to verify the 
eventually gap dependence to the number of neutrons N. The considered Hamilto- 
nian contains the realistic Argonne AV8' NN interaction as well as the Urbana UIX 
TNI; the AV8' finite-size effects as described in Sec. I5.1l and l5.2l were included. 

The FP- AFDMC was used to calculate the energy of low-density PNM by sim- 
ulating several systems with N=12-^18 and N=62-^68 neutrons in a periodic box. 
The considered Hamiltonian contains both AV8' NN and UIX TNI interactions. 
In some case only the NN interaction was employed. All the results obtained with 
FP- AFDMC for several densities are reported in tables 15.121 15.131 15.141 and 15.151 
In each table is reported the energy per neutron evaluated for each considered N. 
In the case of kF=0-4^ fm^^ all the calculations were performed for both AV8' and 
AV8'-hUIX Hamiltonians, while for /f=0.2 fm'^ kF=0.6 fm'^ and kF=0.8 fm^^ 
full calculations were performed only with the full Hamiltonians. 

In order to verify that the system is in a superfluid phase for N=14 and 66 
neutrons a calculation using a trial wave function of the type of normal PNM was 
performed. The energy of the system in the normal phase is compared with the BCS 
value in table [5TT61 At kp=0.2 fm~^ and kp=0.6 fm~^ the superfluid is well favorite 
compared to the normal state, and all the energies obtained with the BCS are lower 
respect to the normal trial wave function. It is interesting to note as using N=66 
at kF=OA im~^ the favorite state is again the BCS, but for N=14 the two states 
are almost degenerate, being the two energies within error bars. Instead, at kp=0.8 
fm~^ the situation is opposite: the BCS state is favorite for N=14, but using N=66 
it seems that the ground-state of the system is the normal phase. 

The energy gap was calculated using a different number of neutrons. In table [5TT71 
we summarize all the results obtained using FP- AFDMC and reported in the above 
tables, and the available CP- AFDMC gap energies taken from Refs. [501 1172j . As 
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0.359(2) 
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0.354(2) 




65 


0.359(2) 
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0.355(3) 
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0.3(3) 


68 


0.361(4) 








0.3(3) 



Table 5.12: FP-AFDMC energy per particle calculated with different number N 
of BCS-paired neutrons with Fermi momentum kF=0.2 fm~^ interacting with 
AV8'+UIX Hamiltonians in a periodic box. For odd N the energy gap is reported, 
and the mean value of A given by N=62-=-68 is indicated. All the energies are ex- 
pressed in MeV. 
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Table 5.13: FP-AFDMC energy per particle calculated with different number N 
of BCS-paired neutrons with Fermi momentum ^^=0.4 fm"-*^ interacting with 
AV8'+UIX and AV8' Hamiltonians in a periodic box. All the energies are expressed 
in MeV. 

it can be seen some finite-size effects are present. With the full AV8'+UIX at each 
considered density the gap energy calculated using N=62^68 is smaller with respect 

to the N=12-^18 case. Also with the AV8' at /i-^-=0.4 fm^^ this behavior is evident. 
This fact could indicate some finite-size effect due to the volume of the simulation 
box. Although the potential energy of the system is obtained by considering the 
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Table 5.14: FP-AFDMC energy per particle calculated with different number N of 
BCS-paired neutrons with Fermi momentum fci?=0.6 fm~^. See table lSTTSl for details. 
All the energies are expressed in MeV. 

effect due to inter-particle distances in the neighbor cells, some quantum correlation 
could still be missing. However, it has to be noted the opposite behavior given by 
CP-AFDMC using AV8' at ^^.=0.6 fm'^ 

The gap calculated with FP-AFDMC using N=12-^18 is in good agreement with 
the CP-AFDMC results, while in the only comparable case of N=62-^68 with AV8' 
at /CF=0.6 fm-i the FP-AFDMC energy gap is significantly smaller. 

The AFDMC algorithm was applied to microscopically study the ^Sq superfluid 
gap in low density neutron matter using realistic Hamiltonian containing a modern 
NN interaction as well as a TNI. The results are summarized in table 15.181 The 
AFDMC allows to solve the ground-state of the Hamiltonian by taking into account 
all the polarization and medium modification effects that all other theories neglected 
or cannot include in the calculation in a realistic way[173j. 
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Table 5.15: FP-AFDMC energy per particle calculated with different number N of 
BCS-paired neutrons with Fermi momentum kF=0-8 fm~^. See table [5T3] for details. 
All the energies are expressed in MeV. 



Hamiltonian 




Ebcs/A En /a 


AVS'+UIX 


0.2 66 


0.355(3) 0.367(2) 


AVS'+UIX 
AVS'+UIX 


0.4 14 
0.4 66 


1.281(3) 1.280(3) 
1.239(2) 1.289(2) 


AVS'+UIX 
AVS'+UIX 
AVS' 


0.6 14 
0.6 66 
0.6 66 


2.584(2) 2.614(4) 
2.579(2) 2.606(4) 
2.593(4) 2.621(4) 


AVS'+UIX 
AVS'+UIX 
AVS' 


0.8 14 
0.8 66 
0.8 66 


4.254(4) 4.269(4) 
4.305(3) 4.277(7) 
4.261(3) 4.218(3) 



Table 5.16: Energy per neutron of the system in both normal {En/ A) and BCS 
{Ebcs/A) phase for several Hamiltonians, densities and N. The estimates of the en- 
ergy with the system in the normal phase were performed with the same calculation 
described in the PNM section. All the energies are expressed in MeV. 
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Table 5.17: Energy gap calculated with FP-AFDMC {App) and with CP-AFDMC 
(Acp) using different Hamiltonians for several densities and considering different 
numbers of neutrons. Some of the CP-AFDMC results are taken from Ref. [50] and 
other from |172j . All the energies are expressed in MeV. 
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A A/Ep 


0.2 


0.3(3) 


0.55 
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1.5(2) 


0.74 
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2.1(2) 


0.47 


0.8 


1.6(4) 


0.20 



Table 5.18: Energy gap calculated with FP-AFDMC with N=62^68 using the 
AV8'-|-UIX Hamiltonian for several densities. All the values are expressed in MeV. 



Chapter 6 



Results: neutron-rich isotopes 



The structure of stable nuclei is generally dominated by shell effects which can be 
nearly completely explained by the single-particle picture of the system. 

However, it has been found experimentally, by means of fragmentation of heavy 
ion beams [mi 11751 11761 1177] . that neutron rich nuclei show features which seem to 
go beyond the standard shell-model. For example, it was observed that the neutron 
drip line for O and F isotopes has a sudden change. This has been interpreted as a 
clear sign that many-body effects in this regime cannot be treated as a perturbation 
of the one-particle picture, but become relevant in determining the overall structure. 
Theoretical discussions of the insurgence of a new magic number at N=16 have re- 
cently been developed by means of mean field calculations [178 j and antisymmetrized 
molecular dynamics [T79]. 

On the other hand the study of neutron-rich isotopes is an interesting investiga- 
tion field because their structure is particularly relevant to constrain properties of 
the crust in neutron stars [180J, but the impossibility to use accurate methods such 
the GFMC to study properties of medium nuclei in addition to the complexity of 
NN and TNI nuclear forces imposes the need to introduce a simpler model rather 
then the full microscopic description of a nucleus with all nucleons as degrees of 
freedom. 

6.1 Modelization of neutron-rich nuclei 

A Green's Function Monte Carlo calculation was performed by mimicking oxygen 
isotopes as neutron drops [381(35]: in that model all neutrons interact with a realistic 
interaction (AV18-I-IL2), with the addition of an artificial external well providing 
the attraction necessary to produce a bound state of neutrons and representing an 
average effect of the protons which are not explicitly included in the calculation. The 
study of oxygen isotopes is then reduced to the study of the N neutrons in the ^+^0 
nucleus. The chosen external well was a Wood-Saxon potential whose parameters 
were determined in order to reproduce the correct separation energy of ^^O and ^^O 
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isotopes |35). 

The CP-AFDMC was applied to study neutron-rich isotopes of oxygen [5T] and 
calcium [52] in a any more simple model, where only external neutrons are included 
in the calculations, while all the core is modeled with a single-particle potential 
acting on valence neutrons. The advantage of this model is that the external well is 
generated in a self-consistent way, and the absence of parameters implies no need 
for any additional calculation. The study of oxygen isotope excited states were also 
performed. 

The considered Hamiltonian is 

^ = - E + E ^-*(^^) + E % + E ^^^^ ' (6-1) 

1=1 i=l i<j i<j<k 

where the NN interaction is the Argonne AV8', the TNI is the Urbana UIX and the 
N are the interacting external neutrons of the i^+A^ q g^j^j ^*'+^Ca isotopes. The one- 
body potential Vext describes the core of the closed-shell nucleus which is replaced 
by a self-consistent potential. 

The trial wave function used in the AFDMC calculation has the same structure 
of that employed in neutron drop calculations. The radial component were obtained 
solving the Hartree-Fock (HE) problem with the Skyrme forces which also provides 
the external potential well that describes the closed core. The yielded radial functions 
are written in the j, rrij base, so the single-particle wave function were obtained by 
coupling the spherical harmonics with the spin using the Clebsh-Gordan coefficients. 
In "O n=18...22 isotopes neutrons fill only the orbitals in the I-D5/2 shell except for 
the isotope ^^O, where it is necessary to put some neutrons also in the ii^ 
order to build the ground state of correct symmetry. In the "Ca ri=42...48 isotopes 
it is assumed that neutrons fill only the orbitals in the 1Fy/2 and shell in order 

to build the ground state of correct symmetry. 

The many-body states are obtained by coupling the single-particle angular mo- 
mentum by constructing eigenstates of total angular momentum J = ji + ■■■ + Jn] 
for states with an even number of neutrons, the ground-state has J = 0, while for 
odd neutron numbers, the ground state has total angular momentum J = 5/2 or 
J = 7/2. These states are in general written in terms of a sum of Slater determinants 
whose coefficients are determined by the symmetry of the state. Each determinant is 
evaluated at the current values of the positions and spin assignments of the nucleons 
in the walker \R,S). 

The choice of the trial wave function proved to be particularly important for 
calcium isotopes calculations. A sufficiently good representation of the ground state 
of ^^Ca can be obtained by building a two-hole state which is complementary to 
^^Ca. For the ^^Ca and ^^Ca completely different trial wave functions are needed; in 
fact the wave function for the first nucleus contains 9 determinants and the second 
35. In principle we could use the ^^Ca wave function as a three-hole state describing 
"^^Ca. However, AFDMC gives in this case an energy a few MeV too high. This is due 
to the approximation used to deal with the sign problem. The two wave functions 
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are degenerate in the eigenspace of but they give two different energy, because 
they have a different nodal surface. 

The same procedure is used for the construction of the oxygen excited states, 
always limiting the choice of the single particle orbitals to the shell lD^i2- 

In these calculations contributions for the 25^/2 orbital were neglected, despite 
one might expect some mixing of this level in the states of isotopes considered. 
However, the fact that AFDMC gives results for the energy differences in very good 
agreement with experimental data indicates that the effects of the mixing are not 
extremely significant. The inclusion of contributions from the 25^/2 orbital also 
gives rise to major technical difficulties. In fact, the component along this orbital 
tends to be projected over the more bound 15i/2 state, giving rise to a non-physical 
density for the external neutrons. This difficulty also prevented from extending this 
calculations to higher mass isotopes to the drip line. 

An important issue in this calculation is the choice of the Skyrme force param- 
eters to be used to generate the effective-HF potentials included in the Hamilto- 
nian for AFDMC calculations. The use of a self-consistent potential well makes the 
single-particle part of the importance function essentially parameter less. In absence 
of an efficient variational procedure, it was important to establish that this choice 
of orbitals could be in general reasonable, if not optimal. 

Skyrme forces cannot be related directly to the realistic NN interaction used in 
the calculation. One must therefore establish a criterion for discriminating among 
different choices. As a first step the variational energy of the closed-shell isotope 
^^O was computed on the correlated Slater determinants obtained using different 
sets of parameters in the Skyrme Hamiltonian, but using the full Hamiltonian with 
potential AV8'-|-UIX. In particular the Skyrme I and Skyrme II[40j forces were tried. 

The expectation value of the Hamiltonian considerably differs in the two cases; 
the variational energy is -29.9(1) MeV using Skyrme I and -60.5(1) MeV for Skyrme 
II, while the experimental value of the energy for the external neutrons is -34.407 
MeV[130j. This is mainly due to the fact that the two sets of Skyrme's parameters 
give quite different effective-HF potentials, as plotted in Fig. 16. 1[ As a criterion 
for determining which shell (and which set of orbitals) to use, it has been chosen 
to minimize the difference between {H) and the available experimental estimates 
of the binding energy of the nucleus. Consequently the Skyrme I force was chosen 
also for calcium isotopes calculation. However one expects that the choice of the 
self-consistent potential only influences the absolute total binding energy of each 
isotope. The difference between energies of different isotope depends mostly on the 
choice of the NN and TNI interactions that include correlations between the external 
neutrons. 

6.2 Oxygen isotopes 

The AFDMC energies obtained for the isotopes series ^^O-^'^OfSlJ are reported in 
table 16.1] compared with the available experimental values p^O] . 
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r [fm] 

Figure 6.1: Effective single-particle potential generated by HE algorithm, with 
Skyrme's parameter of type Force I (full line) and II[40j (dashed line). 



isotope 


EaFDMC Eexp 


E(^«0) 
E(i90) 
E(200) 
E(2iO) 
E(220) 


-18.04(2) -12.188 
-22.4(1) -16.145 
-29.36(6) -23.752 
-33.6(2) -27.558 
-40.48(5) -34.407 



Table 6.1: Ground-state energy calculated with AFDMC. All the energies are express 
in MeV. 

As expected, the computed values are quite different from the experiment, al- 
though the relative discrepancy never exceeds 30%. Of course this is a drawback of 
the use of an external potential for including the effects of the filled core of the nu- 
cleus. In particular the total binding energies are all overestimated. This reflects the 
absence of a correct description of the density of neutrons at the center of the drop, 
which is underestimated, giving rise to an effective potential which is too deep for 
small distances from the center. Moreover core-polarization effects were completely 
neglected. However, most of the information needed to understand the effects of 
NN interaction in the external shell can be obtained looking at energy differences 
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between the isotopes considered. In fact, if the intrashell interaction has a dominant 
effect, the gaps should not depend too much on the quality of the external well 
considered. In Table 16.21 and in Figure 16.21 are reported the energy differences for 
the isotopes considered, compared with the corresponding differences obtained from 
the experimental results. The AFDMC calculations have also been compared with 
the GFMC results reported in the paper by Chang et al. |38j . As it can be seen, in 
this case the agreement between computed and experimental values is excellent. 





EaFDMC Eexp 


E(^yO)-E(^«0) 
E(20O)-E(i8o) 
E(2i0)-E(i80) 
E(220)-E(i80) 


-4.4(1) -3.957 
-11.32(8) -11.564 
-15.5(2) -15.370 
-22.44(7) -22.219 



Table 6.2: Ground-state energy calculated with AFDMC. The differences between 
the studied isotopes are reported. All the energies are express in MeV. 
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Figure 6.2: Outline of differences between energies of the isotope series studied. All 
the energies are expressed in MeV and all values are referred to that ^'^O. Results 
are compared with experimental results [130] and with GFMC calculations of Ref. 



It is possible to notice how the binding energy of O relative to O computed 
by AFDMC is almost coincident with the experimental finding, while there is a 
strong overestimation in the GFMC calculation. This is quite surprising, being ^^O 
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a closed shell nucleus. However, it is possible that in this case the choice of the 
self-consistent HE potential gives a better description of the outer shell than the all 
neutrons picture used in the GEMC calculations. 

In figure 16.31 the AEDMC densities normalized to unity of the external neutrons 
for the isotopes considered are displayed. As it can be seen the neutron densities are 
all quite similar, and very small deviations are present. In fact the density profile 
appears to be larger when there are more neutrons; this can be expected because 
of the repulsive nature of the NN interaction for neutrons. In the figure it is also 
displayed the density of ^^O calculated with Skyrme interaction, to make evident 
the "halo" effect of the external neutrons. 




1- [fmj 



Eigure 6.3: Densities of external neutrons for all isotopes in the ground-state calcu- 
lated with AEDMC and the Skyrme's density of ^^O. 



By constructing the appropriate wave functions it is possible to study the low- 
lying excited states of a given isotope. Excited states were computed for the isotopes 
^'^O, and ^^O. Calculations are limited to states which can be built while still re- 
maining within the I-D5/2 shell. If the considered states are orthogonal to the ground 
state, the projected wave function will be constrained within the same subspace. The 
AEDMC results compared with experimental results of Norum et al. jlSlj for ^^O 
and of Stanoiu et al. |174j for ^^O are reported in table [631 

In Eig. 16.41 it is shown the results for the ^^O isotope, compared with experiment 
and with multideterminant projection Hartree-Eock (MDHF) results by Morrison et 
al.[182j. The excitation energy of the 2"'" state turns out to be in excellent agreement 
with the experimental findings. The situation for the 4"'" state is much less satisfac- 
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state 


A ZT" A 17 

^i^AFDMC ^^exp 


180 

2+ 
4+ 


1.99(9) 2.0(2) 
5.08(2) 3.6(2) 


210 

1/2+ 
3/2+ 


1.24(9) 1.218(4) 
2.11(9) 2.133(5) 



Table 6.3: Excitation energy calculated with AFDMC of the isotopes O and O, 
compared with experimental results (see text). All energies are in MeV. 



tory. There are essentially two possible sources for this discrepancy. In the case of 
the J = 4 state, the large angular momentum might give the spin-orbit component 
of the potential a relevant role. The presence of spin-orbit terms tends to heavily 
modify the nodal structure of the wave function, and this effect can be due to the 
use of the constrained-path approximation. 
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Figure 6.4: Low-lying excitation energies of the ^^O isotope, compared with exper- 
imental values from ref. |181j . All energies are in MeV. Results are compared with 
MDHF results of Ref. jUSj. 



The excitation energies for the ■^^O isotope are shown in Fig. l6.5i We compare our 
AFDMC results with experiments and with shell model calculations by Brown [183] . 
In this case the agreement between AFDMC results and experiments is very good 
for all the states considered. 
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Figure 6.5: Low-lying excitation energies of the ^^O isotope, compared with experi- 
mental values from Ref. [174] , All energies are in MeV. Results are compared with 
shell model results of Ref. |183j . 



6.3 Calcium isotopes 

The constrained-path AFDMC energies obtained for the isotopes series Cs.-'^^ Cq- 
compared with the available experimental values p!30] are reported in table [67 



isotope 


Eafdmc 


E 


E(4^Ca) 


-25.25(7) 


-19.85 


E(43ca) 


-33.9(4) 


-27.78 


E(44ca) 


-46.3(1) 


-38.91 


E(45ca) 


-52.6(4) 


-46.32 


E(46ca) 


-62.9(3) 


-56.72 


E(47ca) 


-70.2(7) 


-63.99 


E(48ca) 


-80.3(8) 


-73.09 



Table 6.4: Ground-state energy calculated with AFDMC. All the energies are express 
in MeV. Experimental values are referred to ground-state energy of ^'^Ca, taken from 
Ref. [T30] . 

As in the case of oxygen isotopes the absolute binding energies are quite different 
from the experimental results, and the total binding energies are all overestimated. 
In table 16.51 and in Fig. 16.61 the energy differences for the isotopes considered are 
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reported, compared with the corresponding differences obtained from the experi- 
mental results. As it can be seen, also in this case the agreement between computed 
and experimental values is excellent. Small deviations are present only for the two 
cases of isotopes '^^Ca and ^^Ca. 





Eafdmc 


E 


E(4^Ca) 


-E(4«Ca) 


-55.1(8) 


-54.09 


E(43ca) 


-E(48ca) 


-46.4(8) 


-46.16 


E(44Ca) 




-34.0(8) 


-35.03 


E(45Ca) 


-E(48Ca) 


-27.7(8) 


-27.62 


E(46Ca) 


-E(48ca) 


-17.4(8) 


-17.22 


E(47Ca)-E(48Ca) 


-10.1(8) 


-9.95 



Table 6.5: Ground-state energy calculated with AFDMC. We report the differences 
between the isotopes we had studied. All the energies are expressed in MeV. 
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Figure 6.6: Outline of differences between energies of the isotope series studied. All 
the energies are expressed in MeV and all values are referred to "^^Ca. 



In this calculation we completely neglected the interaction between shells with 
different orbital angular momentum. This implies that we could verify only the 
accuracy of NN and TNI interaction of nucleons in the F shell. It was shown in 
other calculations that the Urbana IX TNI seems to be unable to correctly describe 
the correlations between nucleons in different shells [63]. However, our results show 
that the energy separation of isotopes with external neutrons filling the shells with 
the same orbital angular momentum L = 3 can be accurately reproduced. 
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In figure [6T71 we report the AFDMC densities normalized to unity of the external 
neutrons for the isotopes considered in this work. In the figure we also display the 
density of ^''Ca calculated with the Skyrme I force. As it can be seen the neutron 
densities are all quite similar, and very small deviations are present. However ex- 
ternal neutrons seem to be very close to the core of ^'^Ca and because of this one 
should expect that the interaction between the core with external neutrons cannot 
be a satisf 
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Figure 6.7: Densities of external neutrons for all isotopes in the ground-state calcu- 
lated with AFDMC and the Skyrme density of "^^Ca. 



Chapter 7 



Conclusions 



In this work I reported and discussed some recent calculations of different nu- 
clear systems. The Auxiliary Field Diffusion Monte Carlo method combined with 
the constrained-path and the fixed-phase approximations, needed to control the 
Fermion sign problem, was described and applied to the study of both finite and 
infinite systems. The fixed-phase approximation, introduced and discussed in this 
thesis, improves all the existing AFDMC results obtained with the previously used 
constrained-path approximation, in particular in dealing with the tensor-r and spin- 
orbit interaction. 

The fixed-phase constraint allows for a better sampling of the random walks 

compared to the constrained path, leading to an overall lowering of the energies. 
Older calculations showed that CP- AFDMC gives good results for nuclei if a simple 
i;4-like Hamiltonian is used, but the imposed constrained cannot correctly build the 
tensor-r correlations. 

Instead, the FP-AFDMC has proved to correctly project out from a simple 
variational wavefunction the tcnsor-r correlations, needed to include the main NN 
physics due to the vr-exchange. It was proved that FP-AFDMC has reached the 
same level of accuracy of GFMC and the most powerful few-body techniques in 
calculating the binding energies of light nuclei both in the closed- and open-shell 
configurations. The improvement due to the choice of the fixed-phase rather then the 
constrained-path is particular evident in the ground-state energy evaluation of the 
oxygen compared with other variational calculations. The FP-AFDMC was applied 
to simulate nuclei up to A=40, but the first estimates of the CPU time needed for 
the calculations indicate that the method will be practical for much heavier systems. 

Another problem was encountered with CP-AFDMC in dealing with the spin- 
orbit interaction. Some discrepancies were observed in the ground-state calculation 
of neutron drops, were a comparison of the spin-orbit splitting with that predicted 
by GFMC points out some difference. The problem with the spin-orbit correlations 
appeared also in pure neutron matter calculations, and was explicitly underlined 
in a comparison between AFDMC and other many-body calculations. Again, the 
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fixed-phase approximation improves the agreement with GFMC in the prediction of 
the spin-orbit spUtting in neutron drops. 

The equation of state of nuclear matter, simulated by means of a periodic system 
with up to 108 nucleons in a periodic box, was calculated with semi-realistic NN 
potential containing tensor and tensor-r forces. The FP-AFDMC results are of 
particular interest and can be used for a benchmark of other many-body theories. 

A comparison with other many-body methods, such as the Brueckner-Hartree- 
Fock and the Fermi Hyper Netted Chain in the Single Operator Chain approxi- 
mation, points out important discrepancies, even in the case with the considered 
Hamiltonian, that is a simplified form containing the main physics of the system 
due to the tensor-r term. However other terms, such the missing NN operators as 
well as the three-nucleon interaction are included in a very approximate way in both 
BHF and FHNC/SOC techniques. 

In particular, the BHF seems to better reach convergence with the addition of 
the three hole-line contribution in dealing with realistic NN interactions; however 
the inclusion of the three-nucleon interaction (that is essential in nucleonic matter 
calculations) is absolutely out of control and some work to estimate its effect should 
be performed. 

The FHNC/SOC approximation is one way to calculate the energy with respect 
to the wave function provided by the Correlated Basis Functions theory. It is a vari- 
ational estimate of the ground-state energy but the uncontrolled introduced approx- 
imations violate the variational principle. This violation is not due to the theory but 
to the method used to calculate the energy of the system. In particular, the inclusion 
of the first class of elementary diagrams in the FHNC/SOC calculations revealed 
that their contributions are particular important and not negligible. The effect of 
other missing terms and of the SOC approximation is however unpredictable. 

One should note that the FHNC/SOC theory was employed to constrain proper- 
ties of the Urbana TNI on the empirical observations on symmetric nuclear matter 
coming from experimental data, so a natural question concerns the accuracy of the 
present Urbana TNI. 

The equation of state of neutron matter has been calculated by means of the 
FP-AFDMC using the semi-realistic Argonnc AV8' NN interaction and the Urbana 
UIX TNI. It has to be noted that although this Hamiltonian cannot reproduce the 
ground-state of light nuclei beyond the s-shell, it is the most used to calculate the 
equation of state of nuclear and neutron matter and then used to constrain properties 
of neutron stars. 

Modern and more complex TNI were proposed by the Illinois group by fitting 
the binding energy of light nuclei, and tested up to the carbon ground-state. The 
AFDMC calculation of neutron drop and neutron matter, using this TNI forces with 
a semi-realistic NN interaction, revealed a dramatic dependence of such forces on the 
density, with the consequence that they cannot be used in a calculation to predict 
neutron star properties. In particular the Svr-exchange contribution in pure neutron 
systems gives rise to a binding of neutrons for densities below the po- 
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Hence the obtained results, even if computed with semi-reahstic interactions, 

strongly indicate that in the region of densities below the bulk nuclear density most 
of the physics described by three-body forces is far to be understood and stress the 
fact that much work must be addressed to the investigation of TNI. In addition, 
many-body forces are likely to play a fundamental role in the understanding of 
nuclear matter properties, particularly at medium and high densities, a region of 
great interest in nuclear astrophysics and heavy ion physics. 

The development of the FP-AFDMC opens several important directions for fu- 
ture investigations. First of all, the addition of the spin-orbit and the TNI forces 
in the Hamiltonian shall permit to explore all the physics of medium-heavy nuclei 
with fully microscopic calculations. 

The first problem is to work out an imaginary time propagator that includes 
the full AV18 NN interaction. Presently, in most of the cases the AV18 is evaluated 
as a perturbation of the AV8'. The assumption is correct by considering that AV8' 
contains most of the contributions of AV18 and was obtained with a reprojection 
by keeping only the more important terms. However the operators appearing in 
the AV18 and not in the AV8' are not included exactly in the GFMC calculations. 
Usually the imaginary-time GFMC evolution is performed with the AV8' and the 
energy is calculated perturbatively with the difference between AV8' and AV18, 
that usually is very small, of the order of fractions of MeV in nuclei. This method 
is used also in the FHNC/SOC calculation where only the first NN correlations 
are included in the variational wave function. However there is not a clear proof 
that such calculation is an upper bound to the true energy, and how good this 
approximation is especially for higher densities. 

The AFDMC was implemented using only a NN propagator in a Vg form with 
the calculation of the AV18 energy in pure neutron matter. Some simple simulations 
suggest that this approximation is valid only in the very low-density limit, but not in 
the general case. For example for 14 neutrons at p <0.12 fm~^ the difference between 
AV8' and AV18 is less than 2 MeV per neutron and is 2.7 MeV and 5.1 McV for 
densities of 0.16 fm~^ and 0.20 fm~^ respectively. On the other hand the calculation 
of energy given by the Vg-part of AV18 or AV8' is completely different also for 
low densities. Considering the Vg propagator only, at /9=0.12 fm~^ the energy of 14 
neutrons with the AV8'+UIX Hamiltonian is 14.12 MeV, while it is 3.60 MeV for 
AVIS+UIX. This means that the extra AV18 terms can not be treated as correction 
to the AV8'. Then some useful indications should concern an accurate study of the 
role of pure NN interaction in neutron and nuclear matter and in nuclei. 

Other considerations should be made about three-body forces. Although at 
present the Illinois forces cannot be used to calculate neutron star properties, they 
could be accurate enough to study properties of heavier nuclei. However, by knowing 
the limitation of actual Hamiltonians, at the present the AFDMC can be used to 
study properties of neutron matter, especially in the low-density regime where the 
superfluid phase, and the TNI effect on the pairing gap, can be accurately investi- 
gated. 
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An interesting open question concerns the range of validity of these effective 
forces depending on the density of the system. It should be pointed out that the 
effect and the role of many-body forces were never explored with other methods. 

The NN interaction is dominant in nuclear matter in the low-density regime but 
the addition of TNI gives a strong contribution to the equation of state, and it is 
necessary to correctly reproduce the equilibrium density that all modern NN poten- 
tials all overestimate. How could then be thought that higher many-body forces are 
negligible, or that their contribution is not important? On the other hand when den- 
sity increases phenomena at subnuclear scale probably starts to play an important 
role, and an effective interaction model should contain more terms that are small at 
lower densities. 

The explicit inclusion of the TNI in the AFDMC algorithm reveals some problem, 
because the Hubbard- Stratonovich transformation naturally applies to Hamiltoni- 
ans where only quadratic operators appear. An attempt to write a propagator by 
introducing more auxiliary field is in progress. The main intrinsic difficulty is that 
an higher number of auxiliary variable to be sampled means higher variance in the 
energy evaluation. 

A very interesting alternative is to consider some fictions excited state of nu- 
clcons (by remembering that TNI are caused by excitations of nuclcons), and a 
corresponding NN fictions-interaction. This new interaction depends on the mass of 
these excitations, and in the infinite limit of such masses this NN potential reduces 
to usual TNI forces. The big improvement is that for intermediate masses it natu- 
rally generates many-body forces at infinite-orders. At the present the AFDMC is 
being used to test the validity of this idea. The starting point is to verify if such new 
interaction exactly gives the same contribution of the Urbana UIX in the low-density 
regime (for example in nuclei) when the excitation masses go to infinity. The next 
step should be to increase the density and study the effect of many-body forces that 
naturally appear with a finite excitation mass. 

Some dreams... 

When the mentioned technical question will be addressed, the AFDMC could be 

used to predict the whole structure of the neutron star, with all the phases of the 
matter at each density. The inclusion of the full NN interaction as well the TNI 
is then of primarily importance, being the starting point to predict properties of 
asymmetric nuclear matter with an unprecedented accuracy, and another exciting 
work to do will be the explicit inclusion of hypcrons in the calculation. For this 
purpose it has to be noted that the existing nuclcon-hypcron and hypcron-hyperon 
interactions have the same structure of NN potentials, and so easy to include in 
the AFDMC algorithm. Other interesting problems (like the possibility of neutron 
near the core of the star to form a lattice) as well the role of many-body forces 
(revealing then a completely new physics in the system) should be investigated. In 
addition, using the AFDMC the full Hamiltonian should be (in principle) computed 
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the ground state of all existing nuclei, in order to have a theoretical prediction of 
the Weizsacker formula directly from NN forces. 

Another very interesting achievement is the simulation of a nuclear system start- 
ing directly from the QCD. This calculation should be a great challenge. At present 
the nuclear physics can not be directly derived from the QCD without the use of 
effective potentials derived from EFT (the lattice QCD is not considered because 
limited to hadron physics). The idea of begin a calculation involving nucleons that 
interact with explicit pions degrees of freedom is very exciting and stimulating for 
the future. In fact the computation of the equation of state of nuclear matter by 
including directly the pions in the calculation is absolutely unimaginable at present. 
Some attempt in this direction is also in progress. 
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